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Abstract 


The  Seismically  Stable  Platform  (SSP)  at  the  Central 
Inertial  Guidance  Facility,  Holloman  AFB,  is  a  dual 
reactionary  mass  isolation  platform  designed  to  remove 
19  ®  g  RMS/Hz  accelerations  and '4  0.02  arcseconds  angular 
position  or  tilt  disturbances  in  all  axes  from  a  test 
environment  intended  for  evaluating  ‘'Third  Generation  Gyro" 
inertial  instruments.  Disturbances  are  removed  by  two 
stages  of  pneumatic  isolators  comprising  a  passive  isolation 
system  and  augmented  by  an  active  control  system  to  cover 
the  operational  bandwidtl  from  IB  ®  Hz  to  100  Hz. 

A  dynamic  model  of  the  SSP  confirmed  the  passive 
vibration  transmissibilities  and  identified  severe 
limitations  on  gain  and  phase  margin  to  the  active  controller 
design.  The  proposed  digital  controllers  identified  the 
SSP  to  be  weakly  controllable.  The  discrete  state 
representation  of  the  SSP  and  control  law  exhibited  numerical 
difficulties  detrimental  to  system  stability. 

This  study  recommends  single  stage  pneumatic  isolation 
or  fluid  isolators  monitored  by  disturbance  parameter 
estimation  schemes. 


ANALYSIS  AND  DESIGN 
OF  A  DIGITAL  CONTROLLER 
FOR  A  SEISMICALLY  STABLE  PLATFOPJVI 


I.  Introduction 


The  purpose  of  this  investigation  is  to  analyze  and 
design  an  active  control  system  for  the  Seismically  Stable 
Platform  (SSP)  at  the  Central  Inertial  Guidance  Test  Facility, 
(CIGTF)  at  Holloman  Air  Force  Base,  New  Mexico. 

An  inertial  test  platform  such  as  the  Seismically 
Stable  Platform  must  have  translation  and  rotational  ground 
vibrations  attenuated  or  isolated  below  the  candidate  test 
instrument's  measurement  sensitivity  thresholds.  An  active 
control ler  should  augment  the  existing  passive  isolation 
system  to  limit  SSP  translations  to  less  than  10  ®g  RMS/Hz 
and  angular  tilts  to  less  than  0.02  arcseconds  over  the 
required  bandwidth  of  10  ®  Hz  to  100  Hz.  (Ref  1,  14). 

To  understand  the  implication  of  these  specifications, 
a  historical  overview  of  seismic  platforms  is  discussed  before 
the  problem  and  approach  are  revealed. 

Historical  Perspective 

Seismic  isolation  platforms  serve  as  a  controlled 
environment  to  test  and  evaluate  inertial  grade  instruments 


needed  in  aerospace  navigation  and  guidance  applications. 

New  missions  have  demanded  greater  accuracy  and  sensitivity. 

To  meet  these  instrument  requirements,  the  Third  Generation 
of  Gyros  (TGG  )  has  emerged  from  instrument  technology  making 
past  test  procedures  and  environments  ineffective.  (Ref  15)  .  To 
judge  properly  instrument  attributes,  the  bias,  drifts,  and 
accuracy  profiles  must  be  distinguishable  and  predictable 
from  test  platform  background  disturbances.  Angular  tilts 
and  specific  force  or  acceleration  disturbances  are  critical 
in  gyroscope  and  accelerometer  sensor  evaluations.  Distur¬ 
bances  may  be  either  measured,  modeled  and  compensated 
in  test  profiles  or  completely  (or  nearly)  removed  or 
isolated  from  testing  parameters.  The  SSP  will  isolate 
disturbances . 

Early  Air  Force  efforts  in  gyro  testing  and  tilt 
stabilization  were  at  the  Frank  J.  Sieler  Research  Laboratory, 
United  States  Air  Force  Academy  using  the  Iso-Pad  isolation 
test  platform  (Ref  21).  Numerous  efforts  (Ref  4,  5,  23) 
were  made  to  control  tilt  or  position  as  well  as  angular 
acceleration  (translational  acceleration,  in  accelerometer 
testing) .  Tilt  control  was  possible  but  translational 
accelerations  were  limited  to  10  ®g  RMS/Hz  above  10  Hz 
only  ( Ref  3 : 18) . 

On  the  basis  of  the  Iso-Pad  control  design  efforts,  a 
seismic  platform  was  needed  v,-ith  a  lower  structure  resonance. 


Engineers  at  FJSRL  and  CIGTF  reviewed  the  dual  reactionary 
mass  concept  and  had  an  isolation  system  constructed  at 
Holloman  A.F.B.  (Ref  24) . 

Background 

The  dual  reactionary  mass  system  was  added  to  an 
existing  seismic  concrete  block  by  the  contractor.  Measure¬ 
ment  Analysis  Corporation  (MAC) .  Concrete  pillars  were 
added  to  the  seismic  block  and  a  steel  box-beam  ring  was 
placed  on  pneumatic  isolators  atop  the  pillars.  On  the 
ring  were  placed  yet  another  set  of  pneumatic  isolators  with 
a  steel  platform  which  would  contain  the  gyro  table  work 
station.  The  pneumatic  isolators  and  steel  structures 
comprise  a  passive  isolation  system,  as  shown  in  Figure  1.1. 
This  system  is  designed  to  respond  as  a  double  low  pass 
filter  above  20  Hz,  limiting  tilt  less  than  +0.02  arcseconds, 
and  accelerations  to  less  than  10  ®g  RMS/Hz.  Recent 
measurements  by  CIGTF  engineers  (Ref  18)  do  not  support  the 
desired  passive  response.  These  discrepancies  are  expanded 
in  Chapter  IV,  SSP  Passive  Response. 

Problem 

The  objective  of  this  thesis  investigation  is  to  analyze 
and  design  an  active  control  system  that  stabilizes  the  CIGTF 
Seismically  Stable  Platform  (SSP)  to  within  ±0.02  arcseconds 
rotation  and  10~®g  RMS/Hz  translation  in  all  axes  over  the 
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Figure  1.1. 


SSP  Passive  Isolation  System 
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bandwidth  of  10  ®  Hz  to  20  Hz.  These  constraints  result 
from  accuracy  and  sensitivities  generic  to  the  TGG  evaluation 
and  are  expanded  in  Chapter  II,  SSP  Functional  Requirements. 

Scope  and  Assumptions 

The  analysis  fully  investigates  the  SSP  dynamics  model 
to  predict  and  to  verify  passive  response  characteristics. 

A  control  law  is  built  on  the  SSP  dynamics  to  meet  the 
specifications  using  Linear  Quadratic  (LQ)  theory.  (Ref  6:415). 

The  Separation  Theorem  or  Principle  (Ref  12)  is  imposed 
to  consider  the  state  estimation  and  controller  design  issues 
separately.  Basically,  the  control  law  considerations  are 
investigated  assuming  the  SSP  sensors  provide  perfect  state 
information.  The  state  representation  looks  at  the  three 
body  dynamics  of  the  seismic  concrete  block  and  two  structural 
members  of  the  SSP  perturbed  about  a  nominal  system  orientation. 
Only  first  order  effects  and  linearized  equations  are  included 
in  the  SSP  perturbation  models.  The  dynamics  development 
clearly  identifies  when  these  approximations  are  made, 
completely  portraying  rotation  and  translation  coupling  modes. 
Finally,  the  study  assumes  that  the  passive  isolation  system 
is  effective  in  removing  rotational  and  translation  disturbances 
from  20  Hz  to  100  Hz.  The  individual  body  structural  reson¬ 
ances  of  the  SSP  are  assumed  to  be  above  100  Hz  (Ref  14:E-3) . 

The  control  law  assumes  influence  over  bandwidth  of  10  ®  Hz  to 
20  Hz. 
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The  FPS  -  foot,  pound,  second  -  units  system  is  used 
throughout  the  study  to  be  consistent  with  other  SSP  reports 
(Ref  14:D-19,  Ref  20).  Mass  units  are  expressed  in  Ib-sec^/ft 
with  the  acceleration  due  to  gravity,  a  constant,  equal  to 
32.2  ft/sec^  (Ref  10) . 

Approach 

The  plan  of  attack  is  divided  in  two  areas  -  dynamics 
model  and  control  law  development. 

The  control  law  specifies  a  Linear  Quadratic  (LQ)  - 
Proportional  Plus  Integral  (PI)  controller  as  the  optimum 
solution  (Ref  6) .  LQ  cost  weighting  emphasizes  system 
states  critical  to  meeting  a  given  performance  specification. 
Proportional  Plus  Integral  control  structures  are  recommended 
typically  (Ref  6:340)  for  disturbance  rejection  and  were 
proposed  by  MAC  for  the  SSP  (Ref  14:D“22).  Discrete  time 
and  actuator  energy  cost  criteria  validate  the  feasibility 
of  the  control  problem  solution. 

Overview 

The  remainder  of  the  thesis  is  divided  into  a  set  of 
controller  requirements,  dynamics  modeling,  passive  response 
and  LQ  design. 

Chapter  II,  SSP  Functional  Requirements,  covers  the 
important  dynamic  model  constraints  and  expected  control 
responses  needed  to  outline  the  active  controller  design. 
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Chapter  III,  SSP  Dynamics,  develops  the  full  six 
degree  of  freedom  equations  for  each  of  the  three  body 
complete  systems. 

Chapter  IV,  SSP  Passive  Response,  examines  various 
stages  of  isolation  for  transmissibility  and  stability 
analysis . 

Chapter  V,  Active  Controller  Design,  develops  the 
truth  model  and  the  reduced  order  state  equations  into  a 
control  law  resulting  in  a  steady  error  within  the  problem 
requirements.  Reduced  order  alternative  models  are  compared 
with  the  system  truth  model . 

Chapter  VI,  Conclusions  and  Recommendations,  summarizes 
the  dynamics  modeling  and  controller  problem  which  leads  to 
recommendations  for  future  work. 
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II. 


SSP  Functional  Requirements 


Introduction 

To  understand  better  the  Seismically  Stable  Platform 
(SSP)  specifications,  a  discussion  of  the  generating  inertial 
instrument's  requirements  and  the  SSP  is  necessary.  In  1975, 
the  Charles  Stark  Draper  Laboratory  commented  on  gyro  testing 
environments  needed  to  support  the  Third  Generation  of  Gyros 
(TGG)  (Ref  15) .  These  considerations  are  the  driving  force 
behind  SSP  development.  With  these  constraints  in  mind,  the 
SSP  system  analysis  considers  tilt  disturbance  constraints 
which  are  important  in  gyro  evaluation,  and  translation 
acceleration  disturbance  constraints  as  the  prime  error  in 
accelerometer  testing. 

SSP  Description 

The  SSP  is  a  dual  reactionary  mass  isolation  platform 
designed  to  have  a  passive  transmissibility  to  seismic 
disturbances  of  -40  dB/decade  above  20  Hz  to  100  Hz 
(Ref  15:5).  The  active  controller  is  concerned  with  the 
transmissibility  from  10  ®  Hz  to  20  Hz. 

The  dual  reactionary  mass  concept  is  explained  best  by 
describing  the  SSP  subsystems  in  Figure  2.1.  The  top  isolation 
level  or  subsystem  is  a  welded  steel  cylinder  with  reinforced 
walls,  suspended  by  its  top  at  four  symmetrical  horizontal 
points  by  steel  boxed  beam  tabs.  This  upper  level  housing 


8 


Upper  Level 


Subsystems 


is  recessed  for  a  gyro  test  table.  Below  the  suspension 
points  are  four  Barry  AL-133-12  pneumatic  isolators,  referred 
to  as  the  "Second  Stage  Isolation"  in  the  study  analysis. 

The  top  steel  structure  and  gyro  table  are  the  "Upper  Level", 
known  as  primary  mass  or  reactionary  mass  in  other  literature 
(Ref  14)  . 

Directly  below  the  second  stage  isolation  is  a  second 
subsystem,  the  "Lower  Level"  which  is  also  identified  in 
other  reports  as  the  secondary  mass  or  intermediate  mass. 

The  lower  level  is  an  octagonal  box  beam  recessed  on  the 
top  and  bottom  to  receive  the  second  stage  and  first  stage 
pneumatic  isolators  as  shown. 

Appendix  A  discusses  the  subsystem  dimensions  and 
physical  attributes.  The  upper  level  mass  is  twice  the 
weight  of  the  lower  level  to  give  a  two  stage  transmissibility 
attenuation  to  disturbances,  with  the  lower  stage  giving 
the  first  low  pass  second  order  response,  and  the  second 
stage  following  a  lower  frequency  cutoff  second  order  response. 
Chapter  IV,  Passive  Response,  fully  explains  the  theoretical 
responses  and  transmissibility  concept. 

The  SSP  is  centered  on  a  concrete  seismic  mass  atop 
concrete  pillars  or  piers.  Future  analysis  will  include 
the  concrete  piers  and  seismic  block  as  simply  the  "pier". 

The  pier  subsystem  is  physically  isolated  below  ground 
level  from  the  laboratory  vibration  environment  by  a  rubber 
grouted  air  gap.  Figure  2.2  shows  the  general  orientation 
of  the  SSP,  the  seismic  pier  and  laboratory.  The  seismic 
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block  subsystem  and  supporting  soil  should  offer  yet  another 
second  order  low  pass  response  concerned  with  frequencies 
higher  than  the  SSP  isolation.  Again,  Chapter  IV 
develops  the  expected  dynamics  further. 

The  supporting  piers,  lower  and  upper  structures  and 
the  pneumatic  isolators  comprise  the  passive  isolation 
system.  How  the  SSP  problem  specifications  and  the  TGG 
specifications  apply  to  the  passive  system  are  examined 
next . 

Third  G eneration  G yro  (TGG)  Specifications 

Quoting  from  Charles  Stark  Draper  laboratory  (CSDL) 
Report,  (Ref  15:1). 

"The  resolution  of  these  new  inertial 
instruments  is  expected  to  be  so  fine 
that  instrumentation  errors  and  uncer¬ 
tainties  must  be  modeled  and  verified 
in  order  to  establish  the  sources  of 
instrument  noise...  the  following 
assumptions  will  be  made  for  performance 
goals  of  the  next  generation  of 
instruments ; 

(  1)  Angular  motion  uncertainty  10  '’meru. 

( 2)  Linear  motion  uncertainty  10  ®g. 

(3)  Measurement  bandwidth  10  z  to 
10  ^Hz." 

Each  CSDL  -  TGG  specification  is  examined  to  determine 
the  implied  constraints  on  the  SSP  dynamics  and  performance 
specifications . 
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Item  (1)  concerns  the  angular  motion  uncertainty  of 
10  ^  meru's  -  milli  earth  rate  units,  or  1.5  X  10~^  arcseconds/ 
second.  This  conversion  is  shown  in  Appendix  B.  The  SSP 
specifications  do  not  include  an  angular  acceleration 
requirement.  For  a  pure  +0.02  arcsecond  tilt  and  using  time 
averaging  (see  Appendix  B) ,  the  angular  acceleration  over 
the  10  ®  Hz  to  100  Hz  bandwidth  results  in  angular  rates 
from  ±  2  X  10  arcseconds/second  to  2  arcseconds/second . 

This  study  was  not  given  a  pure  angular  acceleration  specifica¬ 
tion.  Translation  accelei ations  are  limited  in  all  axes  to 
10“®  g  RMS/Hz. 

Item  (2)  arises  from  a  navigational  error  of  one 
hundred  feet  on  the  Earth's  surface  averaged  over  an  hour 
flight.  Appendix  B  shows  how  the  10  ®  g  measure  is  achieved 
and  demonstrates  the  use  of  average  linear  velocity  and 
acceleration  approximations.  The  SSP  translation  acceleration 
constraint  is  now  10  ®  g's,  but  future  plans  (Ref  20)  include 
a  nano-g  (10  ^  g)  performance  specification. 

Item  (3)  concerns  the  measurement  bandwidth  which  is 
implied  by  the  need  to  observe  the  test  environment  for 
disturbances  lasting  100  second  periods  to  the  full  test 
cycle  of  120  days.  Appendix  B  shows  this  basic  derivation. 
CIGTF  presently  has  a  goal  of  90  days.  The  SSP  bandwidth 
requirement  matches  the  lower  limit  but  exceeds  the  CSDL 
item  (3)  specification  in  order  to  remove  existing  ground 
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motions  abov'e  10  ’  Hz  to  20  Hz;  20  Hz  and  above  is  auqmenteJ 
with  the  passive  isolation  system. 

The  sources  and  characteristics  of  seismic  disturbances 
are  next  discussed  to  understand  the  source  of  accelerat ion 
and  tilt  vibrations. 

Environmental  Disturb aji c e_s 

The  SSP  passive  isolation  system  and  an  active  control 
system  together  must  remove  translational  and  rotational 
disturbances  detrimental  to  inertial  instrument  evaluation. 
Until  the  possible  disturbance  mechanisms  and  physical 
strengths  (magnitude  excursions)  are  understood,  the  control 
issues  in  relation  to  design  specification  are  not  clear. 

Seismic  disturbances  arise  from  distant  earthquakes, 
micro-seismic  waves,  cultural  noise,  low  frequency  tilt, 
acoustical  noise,  test  environment  temperature  changes, 
and  stray  electro-magnetic  fields  (Ref  15:5). 

Earthquakes  can  induce  accelerations  typically  on 
the  order  of  1.6  X  10~®  g  at  0.2  Hz  and  surface  waves  of 
0.4  inches  displacement  at  0.05  Hz  and  100  X  10  ®  g  (Ref  15:5). 
Micro-seismic  waves,  or  microseisms  are  minute  waves  that 
continually  move  through  the  surface  of  the  Earth  and  are 
caused  by  weather  fronts  and  ocean  waves.  The  most  prominent 
and  important  microseisms  are  at  frequencies  0.14  Hz  to 
0.25  Hz  with  accelerations  of  0.8  X  10  ^  g  to  6  X  10  ^  g 
with  amplitudes  from  400  to  1000  micro  inches  (Ref  15:4). 
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Cultural  noise  is  a  manmade  disturbance  caused  from 
vehicle  traffic,  aircraft,  rocket  sled  tracks  and  mechanical 
machinery  such  as  air  compressors  and  air  conditioners. 

Cultural  noise  typically  peaks  at  29  fiz  and  59  Hz  with 
acceleration  of  10  ^  g  (Ref  11:10)  . 

Low-frequency  tilts  may  be  caused  from  building 
temperature  distortions  in  a  predominately  northerly 
direction  (northern  hemisphere)  of  about  20  arcseconds  with 
seasonal  variations  of  90  arcseconds.  This  temperature 
distortion  can  occur  at  rates  of  34  X  10~^  meru  (Ref  15:2). 

Acoustic  wavefronts  within  the  test  area  can  cause 
low-frequency  motions  much  like  barometric  variations 
caused  by  weather  effects  and  sonic  noises  caused  by  jet 
aircraft.  In  a  test  environment,  pressure  and  temperature 
are  usually  monitored  and  controlled;  but  their  effects  are 
very  critical  in  the  operation  of  seismometers  as  well  as 
any  gradients  across  structures  like  the  SSP. 

Electromagnetic  fields  can  cause  small  torque  disturbances. 
These  field  intensities  also  may  induce  very  low  frequency 
signal  noise  effects  into  sensor  measurements  through 
voltages  induced  in  cabling  directly  or  indirectly  through 
ground  loops  in  instrument  amplifiers  or  A/D  converters. 

Seismometer  theory  is  not  addressed  in  this  study 
since  perfect  system  knowledge  is  assumed,  but  the  disturbance 
characterizations  just  mentioned  for  various  phenomena  have 
one  significant  dilemma.  All  the  disturbances  typified  are 
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using  measurement  technology  that  cannot  approach  periods 
of  longer  than  1000  seconds  (10  ^  Hz)  and  the  required 
10  ®  g's  sensitivity  (Ref  8).  Projected  future  test  platform 
bandwidths  for  inertial  testing  require  90  day  (lo”®  Hz) 
test  cycles  or  periods.  Measurement  considerations  have  been 
detailed  in  several  reports  (Ref  14:10-11;  20:15-19;  22), 
but  inertial  testing  needs  to  rely  on  parameter  estimation 
techniques  (Ref  12)  and  the  hope  of  future  seismic  instru¬ 
mentation  technology. 

Another  problem  unique  to  the  SSP  is  vertical  leveling 
of  the  four  support  points  for  the  upper  and  lower  masses 
or  levels.  A  problem  arises  in  getting  all  four  isolators 
adjusted  to  level  each  corner  of  the  SSP  structure  (Ref  2) . 
Each  isolator  has  a  mechanical  limit  switch  preadjusted  to 
keep  isolator  height  constant  under  load  variations.  In 
a  four  point  suspension.  Management  Analysis  Corporation 
(MAC)  feared  that  one  isolator  would  be  above  or  below  the 
nominal  plane  for  a  given  level,  and  the  affected  level 
would  warble  or  tip;  similar  to  the  problem  of  cutting  a 
chair's  four  legs  to  the  same  length.  A  three  point 
suspension  avoids  this  tipping  phenomena,  so  the  South  and 
East  isolator  for  each  stage  have  their  air  supply  regulated 
by  only  one  stage  limit  switch.  The  dynamics  and  controller 
analysis  account  for  this  characteristic.  With  problem 
requirements  outlined,  the  SSP  dynamics  model  is  now 
developed. 
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Ill .  Dynamics  Analysis 


Introduction 

The  SSP  dynamics  are  described  by  linearized  differ¬ 
ential  equations  using  force  and  angular  momentum  relation¬ 
ships.  A  full  six  degree  of  freedom  development  -  three 
translational  and  three  rotational  modes,  considers  perturbed 
states  of  the  upper  and  lower  level,  and  pier  subsystems. 

The  perturbed  states  are  described  about  a  nominal 
coordinate  frame  for  each  of  the  three  bodies  for  the  six 
degrees  of  freedom.  Each  coordinate  frame  uses  a  unique 
subscript  notation  to  identify  clearly  the  SSP  dynamics 
derived  from  a  lumped  spring/damper  network  model.  Once 
the  representative  equations  of  motion  are  obtained, 
translational  and  rotational  cross  coupling  terms  are 
identified.  Before  the  analysis  begins,  a  coordinate 
frame  is  defined  for  each  SSP  structure  along  with 
variables  describing  each  structure's  motion. 

SSP  Coordinate  Frames 

Before  even  a  simple  mass  spring  system  can  be  analyzed, 
all  coordinate  frames,  variables  and  references  must  be  defined. 
With  the  SSP  dynamics,  naming  conventions  are  particularly 
important  because  the  resulting  complex  differential 
equations  exhibit  forms  common  to  much  simpler  mechanical 
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models.  Also,  validity  checks  on  the  equations  are  much 
easier  to  relate  physically  to  the  SSP. 

A  body  coordinate  frame  is  fixed  to  the  upper  level 
(U) /  lower  level  (L)  and  the  pier  (P) .  Figure  3.1  shows 
the  body  frames  and  their  orientation.  A  right-handed 
unit  vector  is  defined  by  Us  =  ui  X  U2  for  the  upper  level 
and  similarly  for  the  lower  and  pier  frames.  Notice  at 
U2  »  I2  »  and  Pa  are  oriented  in  the  navigational  North 
sense,  in  the  nominal  position.  Additionally,  a  nominal 
state  is  assigned  with  Uj  ,  I3  and  p3  oriented  along  the 

/s 

local  vertical  03.  The  e  frame,  the  Earth  frame,  is  taken 
as  inertial. 

Rotation  angles  in  each  frame  are  as  shown,  with 
"^u  '  ‘^u  ®u  rotations  about  the  Earth  frame  ei, 

Bz  and  63  respectively  for  the  upper  level  body.  The  lower 
level  and  pier  rotations  are  defined  similarly.  Appendix 
C  details  the  direction  cosine  transformation  matrices 
necessary  to  express  position  and  velocity  vectors  in  the 
respective  body  frames  to  the  Earth  frame.  The  transformation 
matrices  are  derived  using  small  angle,  first  order 
approximations  as  verified  by  Likins  (Ref  10:102) .  Such 
approximations  assume  the  SSP  oriented  in  a  nominal  state 
and  perturbed  in  small  angular  motions  about  the  respective 
centers  of  mass. 

To  understand  the  importance  of  the  direction  cosine 
transformation  and  to  gain  the  insight  necessary  to  write  the 
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full  six  degree  equations  of  motion,  the  motion  of  a 
single  typical  corner  of  the  SSP  is  analyzed  for  a  small 
perturbation  about  a  nominal  orientation. 

Corner  Perturbation  Derivative 

The  perturbed  state  of  a  body  corner  or  isolator 
suspension  point  is  derived  in  Appendix  D.  Details  of 
derivation  relate  translational  and  rotational  modes  can  be 
more  clearly  understood  from  Figure  D.l.  The  vector  notation 
UE,  which  is  a  position  vector  in  Earth  frame  relating  to 
the  upper  level  (U) ,  easterly  (E)  corner,  is  used  extensively 

O 

to  write  the  equations  of  motion.  Similarly  UE  means  the 
time  derivative,  Earth  frame. 

Now  the  isolator  suspension  points  are  defined  in  matrix 
notation;  the  pneumatic  isolator  characteristics  are 
related  in  vector  components. 

Isolator  Lumped  Parameters 

Before  the  SSP  equations  of  motion  can  be  written, 
the  pneumatic  isolator  is  represented  by  a  directional 
lumped  spring/damper  network  or  conceptually  as  spring/ 
damper  matrices. 

According  to  studies  by  MAC  (Ref  14:C-10),  the  pneumatic 
isolator  has  nearly  the  same  resiliency  in  the  vertical  and 
two  horizontal  axes  for  small  load  deflections.  Essentially, 


20 


the  isolator  is  depicted  as  in  Figure  3.2,  Isolator  Model. 
Notice  the  symbology  used  as  well  as  the  matrix  notation. 
For  example, 


implies  a  matrix  damper  (C)  for  the  first  state  isolation 
(F) .  The  matrix  then  is  expressed  in  Earth  (e)  unit  vectors, 
with  Cpjj  meaning  first  stage  isolation,  X  or  ei  direction. 
Likewise  K  would  be  a  matrix  spring.  The  isolator  charac¬ 
teristics  are  expressed  in  Earth  frame  components  without 
direction  cosine  transformation.  This  assumption  is  made 
because  vertical  characteristics  are  due  mainly  to  the 
pneumatic  isolator  itself,  with  horizontal  characteristics 
due  mainly  to  the  isolator  seal  (elastomer  diaphragm) 
stiffness  (Ref  14:C-10).  Since  the  perturbation  of  the 
isolator  would  not  reorientate  these  characteristics 
appreciably,  no  matrix  transformation  is  done.  Vector 
spring  and  damper  effects  of  the  relative  subsystems' 
displacement  and  velocity  accounts  for  any  compensation 
in  isolator  dynamics  orientation. 
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The  Barry  Controls  AL-133-12  Serva-Levl  pneumatic- 
isolator  characteristics  have  been  studied  by  MAC  (Ref  14:C-53) 
Essentially  the  vertical  characteristics  are  a  function 
of  the  individual  isolator's  load  and  air  supply,  MAC 
determined  the  vertical  natural  frequency  (f^^j)  to  be 
1.8  Hz.  From  the  same  report,  the  horizontal  natural 
frequency  is  4.0  Hz  and  3.8  Hz  for  the  upper  and 

lower  level  isolators,  respectively.  Each  isolator  axis 
represents  a  second  order  system,  comprised  of  a  single 
mass,  spring  and  damper  shown  in  Figure  3.3.  The  differential 
equation,  using  Laplace  operator,  zero  initial  conditions, 
is 


ms^  X(s)  =  -csX(s)  -  kX(s)  -1-  F(s)  (3-2) 

X(s)  =  +  (3-3) 

^  s  +  ^ 
m  m 

The  denominator  of  Eq  3-3  is  in  a  classical  second 
order  form  or  in  generalized  function  of  natural  frequency 
(Wj^)  and  damping  ratio  (c)  • 

s^  +  2i;(jj.,s  +  a)^„  (3-4) 

N  N 
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Figure  3 


m 


X 


m,  mass 

k,  spring  stiffness 
c,  damping  coefficient 
X,  displacement 

F,  applied  force,  positive  x  direction 


.3  Single  Mass,  Spring  and  Damper 
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Equalling  like  order  coefficients  of  Eq  3-3  to  Eq  3-4 
the  following  relations  are  implied: 


c 


c 


c 

m 


4ir^fm 


(3-5) 

u)„  =  2TTf  (3-6) 

N 

f,  frequency  in  hertz 

(  3-7) 


also 


(  3-8) 


k  =  (3-9) 

k  =  (2Tif)2n\  (3-10) 


Using  these  relationships  and  the  MAC  derived  natural 
frequencies,  the  second  stage  isolator  spring  and  damper 
coefficients  are  calculated  using  one  fourth  (four  isolators)  the 
upper  level  mass  ( Ref  14 ;C-53) : 


^VN  =  1-8  Hz  toad  =  nyj/4 


where 


"b 


530.22  lb-sec^ 
ft 


132.56  lb-sec^ 
ft 

(3-11) 
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16955.12  Ib/ft 


SZ 


=  89-95  lb-sec 

ft~ 


=  4.0  Hz  (3-12) 

=  83728.99  Ib/ft 

=  199.89  lb-sec 

Ft~ 


Appendix  E  summarizes  all  the  lumped  spring  and  damper 
parameters  used  in  this  study.  References  are  given  to 
appropriate  pages  in  MAC  study  (Ref  14)  which  give  further 
details  and  derivations. 

Now  the  subsystem  coordinate  frames  and  dynamical 
characteristics  are  described  using  matrix  notation  generic 
to  the  SSP. 


SSP  Equations  of  Motion  -  Six  Degrees  of  Freedom 

The  SSP  has  six  degrees  of  freedom  -  three  translational 
(^if  ^2/  ^3)  and  three  rotational  ,  0,  41)  for  three 
bodies  or  structures  -  the  upper  level  (U )  ,  lower  level  (I.) 
and  pier  (P).  An  equation  is  written  for  each  corner  (such 
as  UE,  UN,  UW  ,  US )  of  each  body.  Then  each  equation  is 
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resolved  into  Earth  components  and  the  translational 
and  rotation  effects  segregated. 

A  lumper  parameter  model  for  the  SSP  is 
shown  in  Figure  3.4.  Various  levels  of  isolation  are 
indicated  for  the  soil  or  ground  (G )  isolation,  first 
stage  pneumatic  isolators  (F)  and  second  stage  pneumatic 
isolators  (S)  using  the  matrix  spring  and  damper  concepts. 

The  spring  and  damper  characteristics  are  expressed  in 
e  unit  vectors,  but  the  spring  and  damper  end  points  are 
vector  points.  For  example,  IE  and  UE  have  and 
connecting  them.  IE,  LN,  LW,  LS  all  are  points  in  the  plane 
representing  the  lower  level  steel -boxed  ring.  Similarly 
the  upper  level  platform  and  pier  seismic  block  are 
described . 

From  basic  spring  and  damper  dynamics  for  one  dimension, 
the  following  equations  result  for  the  spring  (F  )  and 
damper  (F^)  forces. 

Fg  =  kx  (3-13) 

where 

X,  spring  displacement 
k,  spring  constant 

in  vector  case,  using  column  vectors 

Fg  =  [Kg]  i  (3-14) 
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where  F,  x,  K 


■ 

' 

F 

X  = 

X 

X 

9 

F 

y 

y 

F 

z 

2 

. 

>  • 

(3-15) 


^sx 

0 

0 

= 

0 

0 

0 

0 

^SZ 

• 

“  *■ 

^x 

k 

sx 

0 

0 

"y 

= 

0 

^SY 

0 

X  = 

kgyy 

^z 

0 

0 

’^sz 

’^sz^ 

. 

. 

(3-16) 


(3-17) 


The  form  of  Eq  3-14  is  chosen  to  malce  algebraic 
manipulation  easier. 

The  equations  for  the  damper  in  scalar  and  vector  case 


are 


F  =  cx  (3-18) 

c 

where, 

O 

X,  damper  velocity 
c,  damping  coefficient 
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(3-19) 


F 


cs 


[Cg]  X 


where  F,  x,  Cg 


• 

Si. 

o 

■ 

F 

,  X  = 

X 

0  0 

X 

sx 

o 

F„ 

y 

Co  = 

0  c_,,  0 

Y 

s 

SY 

o 

F 

Z 

0  0  Coo 

z 

sz 

^SX 

0 

• 

0 

=s/ 

= 

0 

^SY 

0 

SL 

X  = 

0 

0 

°sz 

o 

. 

■ 

. 

(3-20) 


Again  Eq  3-19  is  written  in  a  form  to  simplify 
equations  later  in  development.  The  displacement  vector 

—  A 

X  and  velocity  vector  x  concept  can  easily  be  written  using 
_  _a- 

UE  and  UE  notation  by  writing  components  as  column  vectors 
in  the  e  frame.  Eqs  D-6,  D-10,  in  Appendix  D,  show  UE 

O 

and  UE  in  vector  notation. 

Extending  the  above  dynamic  force  vectors,  F  =  ma 

can  be  written  for  the  SSP  from  Figure  3.4.  The  SSP  is  a 

constant  mass  system.  Newton's  Second  Lav;  concerns  the 

inertial  center  of  mass  acceleration  of  each  body  as  U, 
as> 

L  and  P  for  the  upper  level,  lower  level  and  pier, 
respectively . 
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Summing  the  forces  on  each  body  due  to  isolator 
spring/damper  equivalents  and  active  controller  actuators, 
results  in  the  following  equations.  For  the  upper  level, 
the  summing  forces  yields. 


m„U 


CglLE 

-  UE) 

0 

+Cg (LN 

-  UN) 

+Cg(I^ 

-  Tm) 

+Cg (LS 

-  ife) 

^UE 

^UN 

+  Kg(LE  -  UE) 
+  Kg(LN  -  UN) 
+  Kg  (LW  -  W) 
+  Kg (LS  -  US) 

F  +  F 
m  US 


(3-21) 


Where  F^^  ,  ,  Fy^ 


and  P, 


US  are  the  four  active 
controller  actuators  between  the  upper  level  and  pier. 
Similarly,  the  lower  level  equations  are  as  follows. 


lOyl.  =  -Cg(LE  - 
-CgiLN  - 

-c^  (rfi  - 
-C3 (is  - 


UE) 

-  Kg(LE  - 

UE) 

-a-. 

-  Kg(LN  - 

UN) 

UN) 

_9_, 

-  Kg (LW  - 

UW) 

UW) 

_a_ 

-  - 

US) 

US) 
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+C„(PE  -  LE)  -  K„(PE  -  LE) 

r  r 

SL.  _  _ 

+Cp(PN  -  LN)  -  Kp(PN  -  LN) 

o  o  _  _ 

+C_(PW  -  LW)  -  K^(PW  -  LW) 
t  r 

+Cp(^  -  LS)  -  Kp{PS  -  LS)  (3-22) 

Likewise,  the  pier  force  equations  are  as  follows: 


-Cp(PE 

-  LE) 

-  Kp(PE 

-  LE) 

O 

a 

-Cp(PN 

-  LN) 

-  Kp(M 

-  LN) 

-Cp 

-  i^) 

-  Kp(PW 

-  LW) 

O 

-  Kp(PS 

-Cp (PS 

-  LS) 

-  LS) 

,-a_ 

_2_ 

- 

-  PE) 

+Cj,(GE 

-  PE) 

,-a_ 

o 

-  ra) 

+C^(GN 

-  PN) 

-  k^(gn 

a 

+C^ (GW 

-  PW) 

-  Kg(^ 

-  PW) 

a 

,  0 

-  PS) 

+C^ (GS 

-  PS) 

-  Kjg(^ 

^PE  “  ^PN  ^PW  ^PS 


(3-23) 


Notice  that  the  -Fp^  ,  “Fp^^  ,  ~Fp^  and  -Fpg  forces 
result  from  the  four  active  controller  actuators  being 
connected  between  the  upper  level  corners  and  the  pier 
foundations . 
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Appendix  F  manipulates  Eqs  3-21,  3-22  and  3-23 
into  differential  equations  for  the  e  components  of  each 
body  acceleration  and  applied  dynamic  forces.  Appendix  F 
yields  the  differential  equations  describing  the  three  degrees 
of  translational  for  the  SSP .  Although  rotational  terms 
are  present,  two  important  points  are  made.  First,  for 
small  angular  movements  the  rotational  terms  certainly 
aid  translational  positions  and  velocities,  as  a  quick 
units  check  verifies.  Secondly,  if  the  SSP  isolator 
placement  is  theoretically  symmetrical,  all  rotation 
contributions  vanish,  which  is  physically  reassuring  showing 
that  the  equations  are  correct.  Now  the  three  degrees  of 
rotational  freedom  are  investigated  to  complete  the  dynamics 
description  of  the  SSP. 

Rotational  equations  of  motions  of  the  three  bodies  - 
upper  level,  lower  level  and  pier,  result  from  summing 
moments  about  each  body  center  of  mass  and  examining 
the  resultant  angular  accelerations.  Expressed  more 
formally,  the  rotational  equivalent  of  Newton's  Second  Law 

_  i  _ 

is  M  =  H  (Ref  10:438),  where  M  is  the  moment  vector 

SL 

about  the  body  center  mass  and  H  is  the  time  derivative  in 
an  inertial  frame  of  the  vector  angular  momentum,  H. 

The  angular  momentum  and  its  derivative  are  expanded 
to  show  the  mathematical  relationship  to  the  body  moments 
of  the  SSP.  Angular  momentum  is  written  as 
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1 


M 


H 


(3-24) 


where , 

H,  angular  momentum 

H,  time  derivative,  inertial  frame 

M,  body  moments 

The  relationship  is  derived  for  the  upper  body  as. 


where , 

ly  ,  moment  of  inertia  matrix 
Uy  ,  angular  velocity  vector 


(3-25) 


(3-26) 


lux  '  ^UY  '  ^UZ  moments  of 

inertia  about  ui  ,  U2  ,  U3  axes. 

000 

+  0U3  (3-27) 

000 

4),  (j),  0  are  tilt  rates  defined  in  Appendix  C. 
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(3-28) 


■ 

■  • 

= 

^ux 

0 

0 

o  ^ 
li'Ui 

0 

^UY 

0 

o  ^ 

(f,U2 

0 

0 

^uz 

o  /v 
0U3 

. 

= 

+ 

0U3 

(3-29) 


The  components  of  the  upper  level  moment  matrix  are 
written  as. 

My  =  My^Ui  +  MyyUa  +  My^Ua  (3-30) 


where, 

^UX  '  ^UY  '  ^UZ  moments  about  ui  ,  ua 

and  U3  axes. 

Equating  components  in  u  frame,  and  transforming  to  e 
frame  using  small  angle,  first  order  approximations  derived 
in  Appendix  C,  the  upper  level  moments  relate  to  angular 
accelerations  as, 


=  lyX^uUi  - 

^^ux'^u 

O  V 

<t>uOu 

'^UY 

- 

(3-31) 

MyyUZ 

oo 

oo 

^^uy'^u 

o  o 

"  ®u^u 

«uz 

Iyjj)]e2 

(3-32) 

MygUa 

v  '' 

oo 

o  o 

*^ux 

-  IyY)]e3 

(3-33) 

35 


Eqs  3-31,  3-32,  and  3-33  are  known  as  Euler's  equations 


(Ref  10:439)  and  may  be  reduced  further  using  small  angle, 


first  order 

approximations  to 

OO 

OO 

( 3-34) 

njY®’ 

\lY  ®  2 

OO 

( 3-35) 

( 3-36) 

U  sing 

similar  transformations  and 

approximations,  the 

lower  level 

is  described  as 

A 

OO 

( 3-37) 

OO 

2 

=  Ily*l®  = 

( 3-38) 

OO 

( 3-39) 

and  for  the 

pier  subsystem, 

OO 

"pxP> 

OO 

( 3-40) 

=  IpY<l>p^2 

OO 

( 3-41) 

■  'pzV’ 

( 3-42) 
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Now  the  body  moments  are  related  to  the  angular 
acceleration  and  moment  inertia.  The  rotation  equations 
are  written  as  vectors  for  each  body  and  respective 
isolator  and  actuator  dynamics.  A  moment  is  easily 
described  in  the  scalar  case  as 

M  =  rF  (3-43) 


where , 

M,  moment 

r,  radial  action  arm 

F,  the  applied  force 

For  the  vector  case,  consider  ax,  y,  z  cartesian 
coordinate  frame,  with  r^^  in  the  X  direction,  F^  applied 
in  positive  Z  and  a  clockwise  moment  about  the  Y  axis. 
In  cross  product  notation  this  reduces  to 

M  =  r  X  F  =  (3-44) 

The  minus  sign  indicates  a  clockwise  moment  for  the 
example.  Now  consider  the  east  corner  of  the  upper  level 
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(UE)  and  the  lumped  parameter  vector  damper  force 
[LE  -  UE]  c  .  Using  the  concept  of  Eq3-19,  the  damper 
force  in  the  e  frame. 

Cg[i^  -  ra]  = 

■‘■^SY^^L  ~  ^U  ^LE®L  “ 

^  \3E®l)^3 

The  moment  arm  of  the  isolator  on  corner  UE  is 

( 3-46) 


86  2 
Cej 

( 3-45) 


Since  the  dynamics  analysis  is  to  be  a  linear  model, 
the  moment  arm  in  Eq  3-46  reduces  to  >  U  se  of  the 

other  terms  would  include  variable  products  and  their 
derivatives  due  to  the  cross  product  operation.  From  a 
physical  standpoint,  under  small  angle  approximations 
and  perturbations  about  a  nominal  body  position,  such 
non-linear  product  terms  would  be  neglectable. 

To  find  the  moment  offered  by  the  UE  damper,  the  cross 
product  is  done  as 
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( 3-47) 


^CS  =  ^ 


Using  Eg  3-47  M^g  is 


M 


CS 


e 1  0  2  e  3 

\jE  ® 

A  B  C 


( 3-48) 


The  indicated  operations  in  Eg  3-48  result  in  the 
components  of  M^g  as 


M 


CS  *  ^Sz\jE^L  ^Sz\jE^  '*■  ^Sz\iE^LE‘*^L  ~  ‘^Sz\jE  \jE ^ 


'‘'^^Sy\jE^L  “  ~  ‘^SY^E^LE®L  “  “^SY^JE  ^E  ^  ^  ®  ^ 


The  first  two  terms  in  the  §2  and  §3  components 
represent  the  translation  cross  coupling  terms  and  the  last 
two  in  each  component  are  the  pure  rotation  contributions. 
Appendix  G  carries  out  the  cross  product  process  similar 
to  Eg  3-48  for  all  SSP  subsystem  corners  and  gives  the 
rotational  differential  eguations.  As  is  done  in  Appendix 
F,  the  eguations  are  separated  into  the  63  ,  §2  and  63  degrees 
of  freedom  with  the  cross  coupling  term  indicated. 
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Summary 


The  differential  equations  developed  in  Chapter  III 
and  tabulated  in  Appendix  F  -  Translation  and  Appendix 
G  -  Rotation  are  the  linearized  differential  equations 
describing  the  upper,  lower  and  pier  level  perturbed 
motions.  These  motions  are  for  a  constant  mass  system 
about  a  nominal  position  in  each  body's  center  of  mass. 

No  acceleration  terms  due  to  gravity  are  present,  because 
the  variance  of  gravity  is  small  and  the  nominal 
orientation  is  an  equilibrium  state  where  all  gravity 
forces,  isolator  forces,  and  actuator  forces  result  in 
the  nominal  position. 

Close  examination  of  the  translational  differential 
equations  in  Appendix  F  clearly  show  the  X  and  Y  directions 
to  have  0  cross  coupling  terms.  The  Z  direction  has  cross 
coupling  terms  from  <p  and  ij^.  These  observations  hold  for 
all  three  levels  and  is  a  physical  check  of  the  described 
motions . 

Similarly,  Appendix  G  shows  the  rotation  differential 
equations  for  each  level  to  have  cross  coupling  present. 

The  ij)  and  <t)  have  Z  cross  coupling  influences  and  0  has  X 
and  Y  cross  coupling  components. 

Chapter  IV  evaluates  the  passive  response  using  Laplace 
operator  solutions  of  the  SSP  differential  equations  in 
the  frequency  domain.  Chapter  V  converts  the  SSP  differential 
equations  into  state  space  notation  for  time  response  analysis 
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using  LQ  synthesis.  The  next  chapter  highlights  the 
SSP  passive  response  to  gain  insight  into  the  actual 
controller  problem  analysis  in  Chapter  V. 

The  dynamics  equations  written  for  the  SSP  give 
insight  to  the  "leveling  problem"  for  adjusting  isolator 
height.  A  nominal  position  would  imply  that  height  of 
the  isolators  would  already  place  a  level  in  equilibrium. 
For  the  study,  a  nominal  position  is  assumed  and  the 
isolators  adjusted  or  leveled.  The  differential  equations 
in  Appendix  F  and  G  could  be  analyzed  in  terms  of  each 
corner  relative  to  the  level  center  of  mass  to  describe 
any  "leveling"  induced  motion. 
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IV.  SSP  Passive  Response 


Introduction 

Before  designing  an  active  controller,  the  SSP  passive 
response  is  analyzed  to  identify  resonant  frequencies  and 
their  effects  on  magnitude  and  phase  responses.  On  the 
basis  of  the  passive  response  characteristics,  the 
systems  poles  and  zeros  then  give  possible  insight  for 
reduced  order  approximations  in  Chapter  V,  SSP  Active 
Controller  Design. 

Simple  mass,  spring  and  damper  networks  are  investigated 
and  their  characteristic  Laplace  solutions  are  related 
to  the  differential  dynamic  equations  resulting  from  Chapter 
III . 

Once  the  vibration  transmissibility  concept  is 
developed,  the  lumped  spring  and  damper  parameters  are 
entered  into  the  transmissibility  expressions  for  each  of 
the  six  degrees  of  freedom  to  obtain  the  system  passive 
response  with  the  controller  actuators  deactivated.  Various 
isolation  levels  are  modeled  to  identify  dominant  SSP 
dynamics  and  possibly  explain  passive  responses  actually 
measured  by  CTGIF. 

Before  the  simple  mechanical  models  are  studied, 
analysis  assumptions  and  methods  are  presented. 
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Analysis  Method  and  Assumptions 

The  passive  analysis  implies  no  controller  actuator 
forces  and  represents  the  SSP  system  open  loop  transfer 
function.  For  the  passive  response,  the  SSP  symmetry 
cancels  any  translation  or  rotation  coupling,  since  all 
isolators  are  equidistant  from  the  body  centers  of  mass. 

In  Chapter  III,  Dynamics  Analysis,  the  pier  center 
of  mass  is  assumed  in  the  plane  of  the  lower  level  isolators. 
The  pier  is  modeled  as  a  plane  having  a  lumped  parameter 
for  inertia  about  the  Z  direction  (ej)  with  equivalent 
spring  and  damper  matrices  representing  soil  dynamics  in 
the  vertical  direction  on  each  corner.  Another  simplifica- 
tion  is  also  made  for  the  moment  of  inertia  in  the  (ez) 
direction  by  using  the  same  lumped  parameter  given  for  the 
X(ei)  direction.  Spring  and  damper  matrices  in  the  two 
horizontal  directions  -  North  and  South,  East  and  West  -  are 
equated  to  the  values  calculated  for  the  X(ei)  direction. 

These  pier  simplifications  basically  consider  the  X  and  Y 
translation  dynamics  as  well  as  the  ijj  and  (J)  rotation  dynamics 
the  same.  In  contrast,  the  pier  simplifications  are  made  by 
assumption  for  analysis  convenience,  while  the  upper  and 
lower  level  dynamics  are  the  same  for  X,  Y  and  ip,  <j>  directions 
by  symmetry.  The  pier  simplifications  are  detailed  in 
Appendix  E  and  do  not  impact  on  active  controller  design 
because  measured  pier  effects  are  above  controller  bandwidth 
at  200  Hz  (Ref  19) .  To  standardize  the  dynamics  qualitatively, 
the  transmissibility  concept  is  defined  next. 


43 


Transmissibility  is  in  general  a  transfer  function  or 
ratio  between  input  and  output  vibrations  or  displacements 
(Ref  17:16).  Expressed  in  the  frequency  domain,  transmissi¬ 
bility  is  written  as  a  transfer  function 

Xo  (jlD) 

T(ja))  =  -  (4-1) 

Xi (jo) 

where,  as  functions  of  frequency 
Xo (jo) ,  output 
Xi (jo) ,  input 

T(jo)  ,  transmissibility  transfer  function 


Transmissibility  logarithmic  magnitude  and  phase 
angle  for  Bode  plot  representation  are  written  as 


LmT (jo) 


20  log 


Xo  (jo) 
xi (jo) 


dB 


(4-2a) 


T(jo) 


tan"^  ^ 

Re[T(jo)  ] 


(4-2b) 


where , 

LmT (jo), 
T(jo)  , 
Im,  Re  , 


log  magnitude,  base  10,  decibels 
phase  angle,  degrees 

imaginary,  real  part  rectangular  components. 
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vibration  transmissibilities  (Ref  17:16)  can  represent 
displacement  responses,  but  also  velocity  and  acceleration. 

By  assuming  T(joj)  to  be  a  linear  transfer  function  and 
sinusoid  velocity  or  acceleration  inputs,  the  displacement 
derivative  responses  are  easily  obtained.  This  assumption 
is  commonly  made  in  seismology  studies  for  low  frequency 
Earth  movements  which  are  naturally  free  of  high  frequency 
components  (Ref  17:16).  For  the  SSP,  the  acceleration 
transmissibility  responses  are  needed  for  the  translation 
passive  responses  and  are  easily  implied  from  displacement 
transfer  functions.  The  SSP  rotational  passive  response 
is  obtained  directly  from  the  angular  position  of  tilt 
transmissibilities.  Once  transmissibility  transfer  functions 
are  computed  for  the  six  SSP  degrees  of  freedom,  probable 
acceleration  (by  approximation)  and  tilt  responses  are 
available . 

Transmissibility  magnitude  Bode  plots  are  also  a 
convenient  representation  of  power  spectral  density  (PSD) . 
Mathematically  the  two  are  related  as 


PSD( jw) 


10 


Xo ( j-) 

Xi  ( j.j) 


2 


20  log 


Xo  (jo)) 
xi ( jw) 


(4-3) 
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The  acceleration  SSP  specification  for  all  axes  is  a 
PSD  requirement  in  terms  of  root  mean  square  (RMS) .  Accelera¬ 
tion  RMS  calculations  are  not  made  because  of  available 
computer  software,  so  only  peak  transraissibility  values  are 
considered . 

Simple  Three  Mass,  Spring  and  Damper  Networks 

A  simple  three  mass  network  shown  in  Figure  4.1  offers 
considerable  insight  into  the  transmissibility  differential 
equation  solution  of  more  complex  networks  such  as  the  SSP. 

From  Figure  4.1,  the  one  degree  of  freedom,  Z  trans¬ 
lational  differential  equations  can  be  written  for  the 
three  masses  using  F  =  ma .  All  accelerations  are  inertial- 
Earth  frame  with  all  masses  constant,  A  displacement  input 
z-  is  assumed  as  shown.  The  equations  of  motion  for  the 
three  masses  are 


OO  0  o 


"*u^u 

■  's'^L  - 

+ 

^s‘^L  - 

(4-4) 

OO 

o 

- 

^S^^L 

- 

o 

+ 

kpCZp  - 

(4-5) 

OO 

IHpZp 

=  -Cp(Zp  - 

o 

- 

kplzp  - 

«g*'g  - 

o 

^p^ 

+ 

- 

(4-6) 
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niy  ,  upper  mass 

m  ,  lower  mass 

Xj 

mp  ,  pier  or  seismic  mass 
kg,  Cg  ,  second  stage  isolation 
kp,  Cp  ,  first  stage  isolation 

k_,  c„  ,  ground  isolation 

Figure  4.1  Simple  Three  Mass  Network 
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Eqs  (4-4),  (4-5)  and  (4-6)  are  rearranged  and  written 
in  Laplace  notation,  assuming  zero  initial  conditions. 


Zy(s)  [n^s 

^  + 

CgS  +  kgS]  =  Zj^(s)[CgS  +  ]<g] 

(  4-7) 

Z^(s)  [m^s 

"  + 

(Cp  +  Cg)s  +  ( kp  +  kg)]  =  ZpS[CpS  +  kp] 

+  Zjj  (  s)  [CgS  +  kg] 

(4-8) 

Zp(s) [mpS 

2  + 

(c^  +  Cp)s  +  ( +  kp)  ]  = 

^(s)[c^ 

,  s  +  ]<^]  +  Zj^(s)[CpS  +  kp)  I 

(4-9) 

Eqs  ( 4-7) ,  ( 4-8)  and  (  4-9)  are  rewritten  with  variables 

representing  the  bracketed  terms  in  preparation  for  algebraic 

manipulations;  these  equations  follow  term  for  term  below; 

Z^(  s)  [a] 

= 

Zj^(s)  [b] 

(4-10) 

ZjA  s)  [e] 

= 

Zp(  s)  [c]  +  2^  (  s)  [d] 

(4-11) 

Zp(s) [f] 

= 

^  {  s)  [g]  +  Z^(  s)  [h] 

(4-12) 

Solving  these  expressions  for  the  transmissability  ^(s)/ 
^(s),  the  transfer  function  is 
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beg 


(4-13) 


^u'=> 


ZqIs) 


aef  -  b^f  -  c^a 


which  is  the  basic  form 


Zy(s)  ^  (Cg  +  kg)  (CpS  +  kp)  (C|gS  +  k^) 

Zg (s)  (m^. . . ) (m^. . . ) (mp. . . ) - (c^s  +  kg) ^ (mp. , . ) 

- (CpS  +  kp) ^ (m^. . .) 


(4-14) 


Eq  4-13  with  the  indicated  operations  carried  out 
to  obtain  a  ratio  of  polynomials,  becomes  a  transfer  function 
of  the  following  form, 

=  As^  +  Bs^  +  Cs  4-  D _ 

Z^(s)  Es®  +  Fs®  +  Gs"  +  Hs®  +  Is^  +  Js  +  K 

(4-15) 

The  polynomial  coefficients  of  Eq  4-15,  A  thru  F  are  tabulated 
in  Appendix  H  in  terms  of  the  variables  for  this  simple 
system  once  the  spring,  damper  and  mass  variables  are 
specified. 

Now  that  a  scalar  Z  transmissibility  transfer  function 
solution  is  available  in  Eq  4-15,  the  translation  equations 
of  Chapter  III  and  Appendix  F  are  reexamined.  Eq  4-4  is 
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repeated  again  in  differential  equation  form  (terms 
collected) . 


Now  consider  the  m^  component  of  03  in  Eq  F-11, 
no  rotation  coupling  terms,  and  actuator  forces  zero. 

00  o  o 

~  ~^^SZ^U  ~  '''  ^^SZ^L  (4-17) 

The  scalar  development  of  three  mass  networks  represented 
in  Eq  4-16  characterizes  the  forms  given  in  Appendix  F. 

Namely,  the  simple  Cg  second  stage  damper  parallels  the  SSP 
second  stage  by  a  factor  of  four,  which  is  physically 
reassuring,  since  the  SSP  second  stage  has  four  isolators. 
Extended  further.  Appendix  H  gives  the  solution  for  three 
stages  of  isolation,  no  rotational  coupling,  no  active 
control  actuator  forces  by  the  simple  substitution: 


^'^SZ 

-> 

^S 

(4-18a) 

^s 

(4-18b) 

^°FZ 

'^F 

(4-18c) 

-> 

kF 

(4-18d) 
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(4-18e) 


4c 


GZ 


4k 


GZ 


(4-18f) 


The  same  parallel  differential  equation  structures 
exist  for  the  direction  in  Eqs  F-5,  F-6,  and  F-7. 
Similarly,  the  ^2  direction  has  the  same  solution  form  by 
using  Eqs  F-8,  F-9  and  F-10. 

For  the  rotational  differential  equations  developed 
in  Appendix  G,  the  parallel  solution  can  be  extended  again. 
Eq  4-16  is  repeated. 

00  O  O 

~°S^U  "  (4-16) 

If  Eq  G-9,  the  component  of  0  is  rewritten  with 
^SX  =  ^SY  isolator  lengths  1^^^  =  l^g  =  1^^  =  1. 

because  of  symmetry,  the  following  rotational  equation 
results . 


00 

'uz°u 


(4-19) 

'""sxi’^u 

Now  a  solution  for  the  3  direction  is  readily  available 
by  simple  substitution. 
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^UZ 

•*> 

mu 

(4-20a) 

(4-20b) 

(4-20c) 

-> 

^F 

(4-20d) 

(4-20e) 

-> 

(4-20f) 

-> 

•^G 

(4-20g) 

Physically  this  follows  from  the  SSP  dynamics  because 
the  Cgj^  ,  ,  CgY  and  vector  springs  and  dampers  act 

as  four  individual  components  about  the  63  axis  for  0  in 
the  upper  level  rotational  dynamics.  Similar  analogies 
can  be  made  for  the  lower  level  and  first  stage  isolation 
effects. 

These  substitutions  do  not  hold  for  and  (}>  .  The 
lyjj  component  of  il;  ,  Eq  G-7  is  rewritten  using  the  same 


assumptions  used 

^uz  ' 

00 

0 

^ux'^u 

t  2kg^l^^ 
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The  ip  substitution  becomes: 


^UX 

-► 

% 

(4-22a) 

-> 

(4-22b) 

2^^821^ 

*^8 

(4-22C) 

2^F2^^ 

^F 

(4-22d) 

2kFzl' 

’^F 

(4-22e) 

2^G2^^ 

-► 

(4-22f) 

2^^021^ 

^G 

(4-22g) 

The  rotational  spring  and  damper  equivalents  are  halved 
because  the  c„„  and  k  dynamics  act  in  pairs  on  opposite 
corners  in  a  see-saw  fashion  about  ei  for  and  §2  for  (}> 
in  the  upper  level.  Similar  comparisons  can  be  made  for  the 
lower  level  and  pier  rotational  equations. 

Deriving  the  Eqs  4-7  through  4-16  is  only  a  vehicle 
for  solving  the  three  mass  isolation  networks.  The  parallel 
made  between  Eqs  4-16  and  4-17  gives  very  realistic  insight 
into  solving  the  @3  equations  found  in  Appendix  F.  When 
the  0,  and  <p  differential  equation  coefficient  substitu¬ 
tions  are  made,  these  scalar  variables  become  rotational 
springs  and  dampers  and  the  Z  notation  of  Eq  4-16  represents 
just  a  scalar  variable,  with  no  translation  implications. 
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SSP  Full  Isolation 


With  the  insight  and  substitutions  gained  from  three 
simple  mass,  spring  and  damper  networks,  the  SSP  full 
isolation  offered  in  the  X(§i),  Y(d2)  and  2(^3)  translation 
directions  as  well  as  the  (J)(d2)  and  0(63)  are 

evaluated.  The  transmissibility  concept  will  provide 
transfer  functions  in  the  required  directions  using  the 
physical  spring  and  damper  values  in  Appendix  e,  the 
polynomial  solution  in  Appendix  H  and  the  substitutions 
just  derived  by  the  simple  mechanical  network  analysis. 

Each  transmissibility  is  defined  in  the  frequency 
domain  from  the  solution  polynomial  and  expressed  in  a  factored 
form  and  a  Bode  plot  from  TOTAL  (Ref  9) .  Actual  polynomial 
coefficients  are  listed  in  Appendix  I.  Table  4.1  summarizes 
break  frequencies,  phase  margins,  gain  margins  and  general 
low  pass  attenuation  offered.  An  expression  for  each 
SSP  directional  transmissibility  is  now  listed  along  with 
an  accompanying  plot  reference. 

Z-Response  (Full  Isolation) 

TZy^^(jw)  =  Zy  ( jw) /Zg  ( ju)) 

=  223.3(10)  +  128.5)  (jo)  +  188.5)^ 

(jo)  +  0.1882  ±  j8,421) 

X  (jo)  +  1.715  ±  j25.37) 

X  (jo)  +  86.47  ±  jl21.4)  (4-23) 

Bode  plot  is  shown  in  Figure  4.2. 
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Table  4.1  Full  Isolation  Frequency  Characteristics  Sununary 
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RESPONSE  (FULL  ISOLATION) 
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FREQUENCY  (HERTZ) 


X-Response  (Full  Isolation) (also  Y-Response) 


TXy/^(ja))  =  Xy(ju>)/X^(jw)  =  Yy(jw)/Y^(ju)) 

=  6Q6.8(ia)  +  199.6)  (ji.0  +  397.9)  (jm 

(jw  +  0.4080  ±  jl8.21) 

X  (jw  +  3.718  ±  j54.94) 

X  (ju)  +  50.17  ±  jl32.5) 

Bode  plot  is  shown  in  Figure  4.3. 


Theta-Response  (Full  Isolation) 

T0u/(,(ju)  =  Qu(  ju)) /0(^(  jw) 

=  (886  .4  (j(.j  +  199.6)  (ju)  +  397.9)  (1 

(jw  +  0.8367  t  jl26.07) 

X  (jw  +  15.65  +  j78.52) 

X  (ju)  +  8.684  +  j81.78) 


Bode  plot  is  shown  in  Figure  4.4. 


PSI-Response  (Full  Isolation) (also  PHI  Response) 


Tlj^y/cijw)  =  (  jU))  /lj)g  (  jw)  =  (  jU>)  /(})q  (  jw) 

=  590.1  (10)  +  128.5)  (iio  +  188.5)^ 

(ju>  +  0.42  ±  jl2.57) 

X  (jo)  +  4.159  +  j39.37) 

X  (ju)  +  42.23  ±  j95.24) 

Bode  plot  is  shown  in  Figure  4.5. 


+  418.9) 


(4-24) 


t  418.9) 


(4-25) 


(4-26) 
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RESPONSE  (FULL  ISOLATION) 


(saay^aoi  uihs  asBHd 

o  o  o  o  o  o 

CD  CO  <o  O  to  CO  (Dm 

m*  ^  I  *-1  ^ 


FREQUENCY  (HERTZ) 


THETfl  -  RESPONSE  (FULL  ISO 


FREQUENCY  I  HERTZ) 


PSl  -  RESPONSE  (FULL  ISOLRTION) 


FREQUENCY  (HERTZ) 


Full  Isolation  Analysis 


The  theoretical  magnitude  and  phase  plots  are  plotted 
using  Eqs  4-2a  and  4-2b  respectively.  Tabulated  values  for 
Table  4.1  are  obtained  from  TOTAL  to  the  indicated  magnitude, 
and  frequency  resolutions.  Analysis  conclusions  are  stated 
now  for  each  direction. 

The  full  Z  isolation  position  resonances  are  theoretically 
at  1.3  Hz  and  4.0  Hz  which  compares  closely  with  the  vertical 
acceleration  resonances  at  1.2  Hz  and  3.5  Hz  measured  by 
FSJRL  engineers  (Ref  19).  A  phase  shift  is  at  2.0  Hz  with 
a  measured  180°  phase  shift  at  about  1.0  Hz  measured  at  CIGTF 
(Ref  19) .  These  good  correlations  give  validity  to  the  SSP 
dynamics  modeling  approach  for  acceleration  transmissibilities . 

No  measurements  are  available  from  CIGTF  for  X  and  Y 
acceleration  transmissibilities.  Figure  4.3  predicts  the 
X  transmissibility  as  well  as  the  Y  transmissibility  due  to 
the  symmetry  of  lumped  spring  and  damper  dynamic  matrices. 
Theoretical  horizontal  position  resonances  are  at  2.9  Hz  and 
8.6  Hz. 

The  0  direction  has  azimuth  resonances  predicted  at  4.1  Hz 
and  12.0  Hz  but  no  measurements  are  available.  Transmissi¬ 
bilities  for  and  cj)  are  the  same  due  to  symmetry  with 
resonances  predicted  at  2.0  Hz  and  12.0  Hz.  No  tilt  trans¬ 
missibility  measurements  are  available  for  the  SSP.  The  0 
direction  has  a  small  gain  margin  and  phase  margin  while  the 
direction  has  gain  and  phase  margins  comparable  to  the  Z  and 
X  directions.  Table  4.1  summarizes  these  points.  No  input 
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signal  magnitudes  are  specified  because  the  transfer 
function  is  calculated  for  a  normalized  response. 

Time  response  characteristics  from  TOTAL  resulting  from  a 
unit  step  input  are  given  in  Table  4.2  for  all  full  isolation 
transmissiblity  functions.  These  time  responses  only  emphasize 
the  lightly  damped  dynamics  of  the  SSP.  Unit  step  responses 
specify  no  input  magnitude  since  they  also  can  be  scaled  or 
normalized  for  a  linear  system 

Simple  Two  Mass,  Spring  and  Damper  Network 

Similar  to  the  three  mass  solutions  just  derived  for 
the  full  isolation  transmissibil ity  transfer  function,  a  two 
mass  system  is  analyzed.  Two  mass,  spring  and  damper  models 
characterize  transmissibilities  result  from  the  three  SSP 
levels  of  isolation  in  active  isolation  combination  two  at 
a  time.  In  actual  operation,  the  SSP  could  have  one  level 
of  pneumatic  inoperative  or  the  SSP  dynamics  could  be  evaluated 
neglecting  the  pier  isolation.  The  possible  scalar  combinations 
for  the  Z  direction  are  shown  in  Figure  4.6.  Again,  the 
whole  driving  force  behind  the  scalar  development  is  to 
characterize  the  differential  equations,  use  the  substitutions 
recognized  in  Eqs  4-18,  4-20,  and  4-22,  and  offer  solutions 
for  the  needed  transmissibil ity  transfer  function.  The 
scalar  Z  is  used  because  it  is  easier  to  visualize  without 
complex  free  body  diagrams  and  it  parallels  easily  into 
the  SSP  Z  direction  solution. 
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Table  4.2  Full  Isolation  Unit  Step  Time  Response  Summary 


Figure  4.6A  SSP  Isola 


Pigure  4.6b  First  Stage/Pi 


6 


upper  mass 


m  ,  lower  mass 

Li 

rUp  ,  pier  or  seismic  mass 
kg,  Cg  ,  second  stage  isolation 
k_,  c_  ,  first  stage  isolation 

r  t 

k^,  c_  ,  ground  isolation 


Figure  4.6C  Second  Stage/Pier  Isolation 


SSP  Isolation 


Figure  4.6A  shows  no  pier  dynamics  or  just  the  SSP 
structure  upper  and  lower  levels.  The  F  =  ma  differential 
equations  for  the  two  masses  are 


oo 

O 

- 

- 

(4-27) 

oo 

=  Cp(Zp 

o 

- 

+  kp(Zp 

o 

- 

-  ^S<^L 

- 

(4-28) 

Rearranging  Eqs 

4-27 

and  4-28 

and  using 

the  Laplace 

operator,  zero  initial  conditions,  and  some  algebraic 
manipulation,  the  transmissability  function  becomes 


Zu(sl 


Zp(s) 


as^  +  bs  +  c _ 

ds"*  +  es^  +  fs^  +  gs  +  h 


(4-29) 


where  the  numerator  and  denominator  polynomial  coefficients 
are 


^  =  ^F^S 
b  =  ‘^F^S  ^S^F 


c  =  kpkg 


(4-30a) 

(4-30b) 

(4-30c) 
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(4-30d) 

(4-30e) 

my(kp  +  kg)  +  CpCg+  nip^kp 

(4-30f) 

b 

(4-30g) 

c 

(4-30h) 

A  solution  similar  to  Appendix  H,  the  full  isolation 
solution  is  now  available.  To  solve  all  four  transmissibility 
functions  -  Z,  X(or  Y)  ,  0  ,  and  i|;(or4))  substitutions  from 
Eqs  4-18,  4-20  and  4-22  are  used. 


First  Stage/Pier  Isolation 

Figure  4.6B  shows  the  second  stage  isolation  inactive. 
The  configurations  have  a  solution  of  the  form  given  in 
Eq  4-29  with  the  polynomials  described  by  Eq  4-30.  If  the 
following  scalar  variable  solutions  are  made,  the 
transmissibility  transfer  function  is  easily  obtained. 


(4-31a) 

(4-31b) 

{4-31c) 

{4-31d) 

(4-31e) 

(4-31f) 
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Now  the  scalar  change  of  variables  are  transformed, 
the  same  substitutions  made  in  Eqs  4-18,  4-20  and  4-22 
are  used. 

Second  Stage/Pier  Isolation 

Figure  4.6C  shows  the  first  stage  isolation  inactive. 
The  following  scalar  variable  substitutions  are  made  in 
the  polynomials  described  in  Eq  4-29. 

+  m  m 

J-i  ir  Lt 

‘'g  *  "f 

°G  °F 

As  is  done  in  the  other  configurations,  the  same 
substitutions  made  in  Eqs  4-18,  4-20  and  4-22  are  used. 
Transmissability  transfer  functions  in  factored  form  are 
obtained  from  TOTAL.  Polynomial  coefficients  are  listed 
in  Appendix  I.  Bode  plot  references  are  indicated, 
frequency  performance  data  and  unit  step  time  response 
information  are  listed  in  Tables  4.3  and  4.4,  respectively. 


(4-32a) 

(4-32b) 

(4-32c) 
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Table  4.3 


Dual  Isolation  Frequency  Characteristics  Summary 

T(ito) 

GM(f ) 

PM(f) 

f  (MP) 

TZy/p(j.)) 

5-25db(  3.0 

Hz) 

5.71°  (2.0  Hz) 

1.3 

4,0 

Hz(  23.57dB) 
Hz(-0  .43dB) 

TZ^/p(ju)) 

39.68db(23.0 

Hz) 

9.42°  (2.6  Hz) 

1.8 

22,0 

Hz(  24.50dB) 
Hz  (-38.93dB) 

TZy^p(ja,) 

38.11db(21.0 

Hz) 

9.42°  (2.6  Hz) 

1.8 

22.0 

Hz(  24.50dB) 
Hz(-38.11dD) 

TXu/p( jw) 

5.29db(  6.5 

Hz) 

5.39° (4 .5  Hz) 

2.9 

8.6 

Hz(  27,99dB) 
Hz(-  0.46dB) 

TXL/p(ja.) 

18.62db(14.0 

Hz) 

8.90°  (5.5  Hz) 

3.8 

16.0 

Hz(  24.63dB) 
Hz(-20.24dB) 

17.46db{14 .0 

Hz) 

8.70°  (5.8  Hz) 

4.0 

20.0 

Hz(  24.69dB) 
Hz  (-22.08dB) 

TOu /p  ) 

6.32db(  9.8 

Hz) 

7.10° (7.0  Hz) 

4.2 

13.0 

Hz(  24.30dB) 
Hz(-  4.56dB) 

T0L/p(jw) 

0.83db{  8.9 

Hz) 

4.20°  (9.5  Hz) 

5.0 

10.0 

Hz(  15.78dB) 
Hz(-  0.71dB) 

T0u/G(ju)) 

2.07db(  8.8 

Hz) 

22.12°(11.0  Hz) 

5.7 

12.0 

Hz(  23.21dB) 
Hz(-  1.92dB) 

Ti^u/p  ( jw) 

6.33db(  4.7 

Hz) 

8.30°  (3.1  Hz) 

2.0 

6.1 

Hz(  24.45dB) 
Hz(-  4.57dB) 

Tij^L/P 

20.97db(13.0 

Hz) 

13.50°(4.0  Hz) 

2.7 

14.0 

Hz(  21.04dB) 
Hz(-22.99dB) 

21.22db(12.0 

Hz) 

13.40°  (3.9  Hz) 

2.7 

13.0 

Hz(  20.99dB) 
Hz(-22.25dB) 

where , 

GM, 

gain  margin 

PM, 

phase  margin 

MP,  magnitude  pea)^ 

AS,  attenuation  slope  approximate  leaving  MP,  towards 
positive  increasing  frequency 


T(jcA)),  transmissability  frequency  function 
(Ref  6:82) 
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Table  4.4  Dual  Isolation  Unit  Step  Time  Response  Summary 


I 


f (t) ,  transmissability  time  function 


Z-Response  (SSP  Isolation) 


TZy/p(ja)) 


=  Zy  (  jw)  /Zp  (  jcj) 

=  1.292  (ju)  +  188.5)^ _ 

(ju  +  0.1883  ±  j8.422)  {joj  +  1.716  ±  j25.37) 


(4-33) 


Bode  plot  is  shown  in  Figure  4.7. 

Z-Response  (First  Stage/Pier  Isolation) 

TZLyp(jw)  =  Zj^  ( jto) /Zq  (juj) 

=  117.3(ioj  +  128.5)  (jo)  +  188.5) 

(ju)  +  0.3391  ±  jll.30)(ju)  +  86.47  ±  jl21.4) 

(4-34) 

Bode  plot  is  shown  in  Figure  '4.8. 
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RESPONSE  (SSP  ISOLATION) 
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FREQUENCY  IHERTZ) 


RESPONSE  (FIRST  STftOE/PIER  ISOLRTION) 


3 


FREQUENCY  (HERTZ) 


Z-Response  (Second  Stage/Pier  Isolation) 


Zy(ju))  /  Zg(ju)) 


=  111.2(i(x)  +  128.5)  (1u)  +  188.5) 

(juj  +  0.3392  ±  jll,30)  (jw  +  82  + 


Bode  plot  is  shown  in  Figure  4.9. 

X-Response  (SSP  Isolation) 

TXy/p(ja))  =  Xy(  jio) /Xp(ja))  =  ( j  w) /Yp  ( j  w) 

«  6.061  (10)  +  397.9)  (ju)  +  418.9) 

(jo)  +  0.4079  ±  jl8.23)(ja)  +  3.717 


Bode  plot  shown  in  Figure  4.10 


jll9.8) 

(4-35) 


±  j55.00) 

(4-36) 
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RESPONSE  (SECOND  STfiGE/PIER  ISOLATION) 


FREQUENCY  I  HERTZ) 
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Figure  4.10  X-Response  (SSP  Isolation) 


TXj^(  joj) 


!First  Stage/Pier  Isolation) 


=  X^(  joj) /Xp(  j..,)  =  Yj^(jw) /Yp(ja)) 

=  143.4(ia)  +  199.6)  (ioj  +  397.9) _ 

(jo)  +  0.1763  ±  j23.81)(ju  +  50.17  ± 


Bode  plot  shown  in  Figure  4.11. 

X-Response  (Second  Stage/Pier  Isolation) 

TXu/G^jw)  =  Xy(  jo)) /Xg(  jw)  =  Yy  ( jcj) /Yg(  jco) 

=  143.2(ia)  +  199.6)  (io)  +  418.9) 

(jo)  +  0.7543  +  j25.08)  (ju  +  47.55 


Bode  plot  is  shown  in  Figure  4.12. 


jl32.5) 

(4-37) 


+  jl29.4) 

(4-38) 
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RESPONSE  (FIRST  STflOE/PIER  ISOLATION) 
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FREQUENCY  (HERTZ) 


Theta-Response 


T0y/p(jU3) 


(SSP  Isolation) 

=  0y( jw) /Op ( j j) 

=  27.66  (io  +  397.9)  (jo)  +  418.9) _ 

(jo)  +  0.8363  ±  j26.10)  (jco  +  8.273  ±  j81.81) 


(4-39) 


Bode  plot  is  shown  in  Figure  4.13. 

Theta-Response  (First  Stage/Pier  Isolation) 

T8j^/p(ja))  =  0j^(ju)/0p(jj) 

=  96.53(ja)  +  199.6)  (iuj  +  397.9) 

(ju)  +  1.508  ±  j34.50)(ja)  +  16.06  +  j78.57) 

(4-40) 


Bode  plot  is  shown  in  Figure  4.14. 
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THETfl  -  RESPONSE  ( SSP  ISOLRTION) 
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FREQUENCY  (HERTZ) 


THETfl  -  RESPONSE  (FIRST  STflOE/PIER  ISOLATION) 
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FREQUENCY  (HERTZ) 


Theta-Response  (Second  Stage/Pier  Isolation) 


0y  ( jw) 


=  97.34  (ia)  +  199.6)  (io)  +  418.9) 

(jo)  +  1.532  ±  j35.71)  (jo)  +  15.92  ± 


Bode  plot  is  shown  in  Figure  4.15. 

PSI-Response  (SSP  Isolation) 

TV^U/p(jw)  =  I/Jy  ( jw) jw)  =  (})y  ( jw) /(fip  ( jto) 

=  6.99  (ico  +  188.5)  ^ _ 

(jw  +  0.4201  ±  jl2.58)  (ju)  +  4.160 


Bode  plot  is  shown  in  Figure  4.16. 


j78.23) 

(4-41) 


±  j39.38) 

(4-42) 
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Figure  4.15  Theta-Response  (Second  Stage/Pier  Isolation) 


PSI  -  RESPONSE  (SSP  ISOLATION) 
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FREQUENCY  (HERTZ) 


PSI-Response  (First  Staqe/Pier  Isolation) 

( jw) /(^p  ( joj) 

=  131.1(10)  +  128.5)  (j  +  188.5) 

(jo)  +  0.7764  +  jl7.09)  (jo)  +  42.23  ± 


Bode  plot  is  shown  in  Figure  4.17. 


PSI-Response  (Second  Stage/Pier  Isolation 

=  i|;y(jo)) /t()g(jo))  =  (})y  ( ju) ( jo)) 

=  125.2(ja)  +  128.5)  (jo)  +  188.5) 

(jo)  +  0.7496  ±  jl6.79)(ju)  +  41.77 


Bode  plot  is  shown  in  Figure  4.18. 


j95.24) 

(4-43) 


±  j94.82) 

(4-44) 
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PSi  -  RESPONSE  (FIRST  STfiGE/PIER  ISOLRTION) 


FREQUENCY  (HERTZ) 


PSI  -  RESPONSE  (SECOND  STfiOE/PIER  ISOLATION) 
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FREQUENCY  (HERTZ) 


Dual  Isolation  Analysis 


Two  important  results  are  gained  from  the  analysis  of 
the  dual  isolation  configurations.  First,  the  SSP  isolation 
dynamics  do  dominate  the  passive  transmissibility  responses. 
Secondly,  the  small  gain  and  phase  margins  are  directly 
attributable  to  the  SSP  isolator  dynamics. 

Table  4.3  summarizes  the  frequency  characteristics  for 
the  dual  isolation  response.  The  Z  direction  is  a  good  example 
of  the  SSP  isolation  dynamics  role.  The  SSP  isolation 
dynamics  in  the  Z  direction  compare  directly  to  the  theoretical 
first  and  second  resonant  peaks  discussed  for  the  Z  full 
isolation,  clearly  showing  the  dominance  of  the  SSP  isolation 
over  the  pier  dynamics.  When  the  first  stage  or  the  second 
stage  isolation  are  considered  with  the  pier,  the  pneumatic 
isolator  dynamics  dominate  and  compare  exactly  to  the  isolator 
vertical  design  frequency  in  Appendix  E.  The  second  resonance 
of  either  isolation  stage  and  the  pier  dynamics  is  at  22  Hz. 

The  SSP  isolation  configuration  in  the  Z  direction  has 
the  same  small  gain  and  phase  margin  as  did  the  full  Z 
isolation  case.  Only  when  one  level  of  pneumatic  isolation 
is  considered  does  the  gain  and  phase  margin  increase  as 
shown  in  Table  4.3. 

The  same  SSP  isolation  dynamics  analysis  and  results 
can  be  extended  for  the  X  (or  Y)  ,  0  and  iJ)(or(j))  directions. 

A  summary  of  unit  step  responses  from  TOTAL  for  the  dual 
isolation  configurations  are  presented  in  Table  4.4.  Clearly, 
the  long  settling  times  and  peak  overshoots  are  identified 
with  the  SSP  isolation  dynamics. 
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Simple  Single  Mass,  Spring  and  Damper  Networks 


A  single  mass,  spring  and  damper  single  level  isolation 
system  is  shown  in  Figure  4.19A.  Figures  4.19B  through 
4.19D  show  the  various  SSP  single  level  isolation  config¬ 
urations  possible  using  the  Z  scalar  derivations  used  in 
the  other  simple  three  and  two  mass  mechanical  networks. 

The  position  transmissibility  for  Figure  4.19A  is  written 
in  Laplace  operator  notation,  zero  initial  conditions: 


Z2  (s) 


cs  +  k 


Zi  (s) 


ms^  +  cs  -T-  k 


(4-45) 


Or  in  the  frequency  domain,  the  transmissibility 
is  written  as 

Z2(jw)  =  cjgj  +  k _  (4-46) 

Zi  ( jw)  mjo)^  +  c j  +  k 


The  transmissibilities  for  the  various  single  levels 
are  solved  by  substituting  the  active  SSP  isolation  level 
spring  and  damper  for  k  and  c  respectively.  The  mass,  m, 
simply  becomes  the  total  mass  supported  by  the  isolation 
level.  These  substitutions  are  summarized  for  the  indicated 
figures  as  follows: 
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Figure  4.19A  Simple  Mass 


Figure  4.19B 
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k 


(4-47a) 

(4-47b) 

(4-47C) 


Figure  4.19C 
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k 


(4-48a) 

(4-48b) 

(4-48c) 


Figure  4.19D  m  +  m  +  m 

U  Li  i 
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k 


(4-49a) 

(4-49b) 

(4-49c) 


Once  the  scalar  substitutions  have  been  made,  appro¬ 
priate  parallel  forms  are  identified  as  s  done  for  the 
other  networks  using  Eqs  4-18,  4-20  and  4-22. 

The  three  transmissibilities  just  solved  are  now 
expressed  in  factored  form  for  the  six  degrees  of  freedom 
with  Bode  plot  references  given.  Tables  4.5  and  4.6 
summarize  key  frequency  and  unit  time  response  performance 
values,  respectively. 
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Table  4.5 


Single  Isolation  Frequency  Characteristics  Sununary 


T  ( jo))  Resonant  Frequency  Peak  Magnitude 


TZu/l^3w) 

1.8  Hz 

24 .45dB 

TZ^/p(j(.) 

1.8  Hz 

24.45dB 

TZp/G(ja3) 

19.0  Hz 

3  .OldB 

tVl'3"' 

4.0  Hz 

24 .45dB 

3.8  Hz 

24 .45dB 

19.0  Hz 

5.42dB 

™u/L<3“) 

5.7  Hz 

21 .39dB 

5.5  Hz 

21 .25dB 

T0p/G(ju3) 

20.0  Hz 

8 .79dB 

2.7  Hz 

20.78dB 

T4'L/p(j<^) 

2.7  Hz 

20.79dB 

T4'p/g( 

15.0  Hz 

4 .37dB 

where , 

T( juj)  , 

transmissibility  frequency 

function 

(Ref  6:82) 


94 


Table  4.6  Single  Isolation  Unit  Step  Time  Response  Summary 
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f(t),  transmissibility  time  function 


Z-Response 


(Second  Stage  Isolation) 


”  Zy  ( jo)) /2;^(  jj) 

=  0.6786(iu)  +  188.5) 

(jo)  +  0.3393  ±  jll.30) 


(4-50) 


Bode  plot  is  shown  in  Figure  4.20. 


2-Response  (First  Stage  Isolation^^ 


TZl/p(3w)  =  Z^( jw) /Zp( j j) 

=  0.6786(10)  +  188.5)  __ 

(jo)  +  0.3393  ±  jll.30) 


(4-51) 


Bode  plot  is  shown  in  Figure  4.21. 
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RESPONSE  (FIRST  STflOE  ISOLATION) 
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FREQUENCY  (HERTZ) 


Z-Response  (Pier  Isolation) 


TZpyG(j(^)  =  2p  ( jo)) /Z^  ( jw) 

=  150.Q(j.jj  +  128.5) 

(jw  +  74.99  ±  jll6.8) 

(4-52) 

Bode  plot  is  shown  in  Figure  4.22. 

X-Response  (Second  Stage  Isolation) 

TXy^j^(ju)  =  XG(jaj)/Xj^(ja))  =  I  jw) /Y^  ( jw) 

=  1.508(iu)  +  418.9) 

(ju)  +  0.7540  ±  j25.12) 


Bode  plot  is  shown  in  Figure  4.23. 


(4-53) 
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FREQUENCY  (HERTZ) 


RESPONSE  (SECOND  STfiCE  ISOLRTION) 
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FREQUENCY  (HERTZ) 


X-Response  (First  Stage  Isolation) 


TXj^/p(ju))  =  Xj^  ( jo)) /Xp  ( jU3)  -  ( juj) /Yp  ( j'o) 


=  1.433  (10)  +  397.9) 

(jo)  +  0.7163  ±  j23.87) 


(4-54) 


Bode  plot  is  shown  in  Figure  4.24. 


X-Response  (Pier  Isolation) 


TXp^^(jo3)  =  Xp(jaj)/XQ(jo))  =  Yp(jw)/Yg(ju)) 

=  86.87  (10)  +  199.6) 

(jo)  +  43.44  ±  jl24.3) 


(4-55) 


Bode  plot  is  shown  in  Figure  4.25. 
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RESPONSE  (FIRST  STAGE  ISOLATION) 


FREQUENCY  (HERTZ) 


RESPONSE  (PIER  ISOLATION) 


Figure  4.25  X-Response  (Pier  Isolation) 


PSI-Response  (Second  Stage  Isolation) 


TiJ^u/l^^w)  =  ijjy  ( jw) (  j  _')  =  (})jj  (  jw) /4)j^  (  jtj) 

=  1 .50  (  jo)  +  188.5) _ 

(jw  +  0.7499  ±  jl6.80) 


(4-56) 


Bode  plot  is  shown  in  Figure  4.26. 

PSI-Response  (First  Stage  Isolation) 

Ttl;L/p(j(o)  =  ij;^  ( joj) /tl^p  ( jw)  =  ( jw) /(fp  ( jco) 

=  1.554  (jto  +  188.5) 

(jo)  +  0.7768  ±  jl7.10) 


Bode  plot  is  shown  in  Figure  4.27. 


(4-57) 
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PSI  -  RESPONSE  (SECOND  STAGE  ISOLATION) 
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PSI  -  RESPONSE  (PIER  ISOLI^TION) 
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THETfl  -  RESPONSE  (SECOND  STfiDE  ISOLATION) 
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FREQUENCY  (HERTZ) 


Theta-Response  (First  Stage  Isolation) 


T0j^^p(ja))  =  0j^(  juj) /0p  {  jo)) 

=  3.012(1ui  379.9) 

(joi  +  1.506  +  j34.59) 

Bode  plot  is  shown  in  Figure  4.30. 

Theta-Response  (Pier  Isolation) 

T0pyG<jw)  =  0p  ( jw) /0G(ju)) 

=  31.32(jaj  +  199.6) 

(ju)  +  15.66  +  j77.51) 

Bode  plot  is  shown  in  Figure  4.31. 


(4-60) 


(4-61) 
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THETfl  -  RESPONSE  (FIRST  STfiOE  ISOLATION) 
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FREQUENCY  (HERTZ) 


THETfl  -  RESPONSE  (PIER  ISOLATION) 


FREQUENCY  (HERTZ) 


Single  Isolation  Analysis 


Analysis  of  the  full  and  dual  isolation  schemes  clearly 
demonstrates  the  dominance  of  SSP  isolation  dynamics.  The 
single  level  isolation  configurations  verify  isolator  design 
considerations  for  translation  resonances  and  yield  some  new 
insight  into  rotational  isolator  effects. 

Reviewing  the  contents  of  Table  4.5,  the  Z  direction 
resonant  frequencies  for  the  second  and  first  stage  isolation 
are  at  1.8  Hz;  the  intended  vertical  resonance  for  each  isolator 
calculated  in  Appendix  E.  Similarly,  the  X  (or  Y)  direction 
has  resonant  frequencies  of  4.0  Hz  and  3.8  Hz  for  second 
and  first  stage  horizontal  single  isolator  resonances  as 
planned  in  Appendix  E. 

An  interesting  result  of  the  rotational  transmissibilities 
is  the  rotational  isolator  dynamics.  The  ©  resonance  is  at 
5.7  Hz  and  5.5  Hz  for  the  second  and  first  stage  resonances. 

A  difference  in  resonance  frequency  for  the  two  levels  is  due 
to  loading  effects  on  isolator  stiffness  as  is  presented  in 
the  X  direction  between  the  second  and  first  stage  isolator 
horizontal  resonances.  The  ’^(or(p)  direction  has  the  same 
resonant  frequency  of  2.7  Hz  for  both  stages  of  isolation 
which  follows  the  same  vertical  frequency  resonance  for  both 
levels  since  the  intended  isolator  vertical  frequencies  are 
designed  the  same.  In  all,  the  translational  and  rotational 
transmissibilities,  the  resonant  magnitudes  are  nearly  five  times 
greater  for  the  second  or  first  stage  isolation  dynamics 
compared  to  the  pier  effects. 


114 


In  general,  the  pier  dynamics  are  prominent  at  20  Hz, 
the  upper  limit  of  the  active  controller  bandwidth.  The 
unit  step  time  responses  in  Table  4.6  show  the  characteristic 
slow  time  responses  introduced  by  either  the  second  or  first 
stage  isolation  dynamics. 

Transmissibility  Pole  and  Zero  Analysis 

The  transmissibility  transfer  functions  for  the  various 
levels  of  isolation  are  analyzed  for  each  of  the  six  degrees 
of  freedom.  For  any  direction,  the  full  isolation  case 
represents  the  complete  set  of  poles  and  zeros.  By 
eliminating  an  isolation  level,  the  effects  on  the  overall 
transmissiblity  can  be  verified. 

The  upper  level  dynamics  in  all  directions  predominate 
the  system  performance.  The  lower  level  poles  are  possible 
sources  of  instability  as  is  discussed  in  Chapter  V,  with 
the  effects  on  the  root  locus.  The  pier  dynamics  do  come 
theoretically  at  the  edge  of  the  active  controller  bandwidth 
at  20  Hz,  but  in  actual  measurements,  the  pier  effects  are 
well  above  200  Hz  (Ref  18) .  Since  the  poles  and  zeros  of 
the  various  isolation  levels  do  parallel  those  of  the  full 
system,  the  solutions  are  verifiable.  The  slight  differences 
in  pole  and  zero  location  are  due  to  scaling  caused  by 
different  configurations  having  different  applied  loads. 

The  rotational  transfer  functions  have  an  interesting 
result.  The  0  transmissibility  zeros  are  the  same  as  those 
of  the  X  or  Y  transmissibility  zeros.  Similarly,  the 
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and  4)  transmissibility  zeros  are  identical  to  the  Z  trans- 
missibility  zeros  because  the  same  spring  and  damper  dynamic 
matrices  are  used  in  the  G  and  X(or  Y)  or  i|)(or(j))  and  X(or  Y) 
transmissibilities . 

The  Z  transmissibility  and  X(or  Y)  transmissibility 
pole/zero  analysis  summaries  are  in  Table  4.7  and  4.8, 
respectively.  For  the  0  and  il)(or(j>)  rotational  transmissibility 
pole/zero  summaries  see  Table  4.9  and  4.10,  respectively. 
Several  analytical  checks  are  possible  using  the  results  of 
these  tables.  First,  the  full  isolation  in  any  given 
direction  is  equal  to  the  combination  of  all  three  single 
isolation  levels.  Secondly,  the  full  isolation  in  any  given 
direction  is  equal  to  two  levels  of  isolation  plus  a  single 
level  of  isolation  in  that  direction.  Care  must  be  taken  in 
these  calculations  because  the  resonant  frequencies  (see 
Eq  3-8)  are  a  function  of  the  isolator  applied  load.  If 
single  levels  are  to  be  added  for  the  full  isolation  verifica¬ 
tion,  each  single  level  must  consider  the  load  applied  to  it. 
Analyzing  Figure  4.1  to  emphasize  the  loading  effect,  the 
second  stage  isolation  load  is  m^  while  the  first  stage  load 
is  mj^  (not  m^  +  nij^)  and  the  pier  stage  load  is  mp  (not 
nip  +  my  +  mj^)  .  Each  pole  and  zero  for  the  single  stage 
correspond  exactly  to  the  full  isolation  poles  and  zeros 
when  properly  loaded.  This  verification  is  reassuring  for 
the  transmissiblities  methods  used  and  the  physical 
intuition  of  the  SSP.  This  verification  was  used  successfully 
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-Transmissibility  Transfer  Function  Zeros  and  Poles 


T(jw)  -  transmissibility  transfer  frequency  function 


Table  4.10  PSI  or  PHI  Transmissibility  Transfer  Function  Zeros  and  Poles 


transmissLbility  tra: 


during  the  early  phases  of  the  study  to  check  the  various 
transmissibility  solutions. 

Besides  an  interesting  analytic  check,  the  poles  and 
zeros  typify  the  resonant  frequencies  and  responses  for  a 
given  isolation  configuration.  These  level  traits  are  used 
extensively  in  Chapter  V  to  build  active  controller  designs. 

Summary 

A  thorough  passive  response  study  for  each  level  is 
not  done  for  two  reasons.  First,  the  pneumatic  isolator 
and  seismic  pier  soil  characteristics  are  not  known  precisely ; 
so  approximating  the  SSP  performance  on  the  basis  of 
equivalent  lumped  parameter  springs  and  dampers  would  only 
be  a  qualitative  treatment.  Secondly,  the  theoretical 
SSP  frequency  responses  can  give  only  general  insight  into 
problems  faced  by  the  active  digital  controller.  Chapter 
V  does  expand  on  the  directional  transmissibility  poles  and 
zeros  using  root-locus  concepts  (Ref  4:203).  Other  important 
points  gained  from  the  passive  frequency  responses  are  now 
summarized . 

The  pneumatic  isolators  do  have  individual  resonant 
frequencies  for  the  upper  and  lower  levels  taken  separately 
(single  stage  isolation)  of  1.8  Hz  for  the  vertical  direction 
and  3.8  Hz  and  4.0  Hz  for  the  horizontal  motion  (single 
stage  isolation)  for  the  first  and  second  stages,  respectively. 
The  full  isolation  theoretically  results  in  the  first  cutoff 
at  1.3  Hz  and  second  at  4.0  Hz.  Actual  SSP  measurement  at 
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Holloman  by  FJSRL  engineers  calculated  (Ref  19)  vertical 
frequency  breakpoints  at  1.2  Hz  and  3.5  Hz. 


Rotational  isolator  dynamics  have  not  been  specified 
by  MAC,  but  theoretical  0  resonant  frequencies  of  5.7  Hz 
and  5.5  Hz  for  second  and  first  stage  isolation  and  2.7  Hz 
for  ip  transmissibilities  (both  stages)  are  predicted. 

In  general,  the  poles  and  zeros  do  analytically  follow 
from  one  (direction)  isolation  configuration  to  another  in 
the  complex  frequency  plane  characteristics,  when  the 
appropriate  load  scaling  effects  are  considered.  The 
transmissibility  solutions  are  mathematically  and  dynamically 
correct. 

Eq  4-35  is  repeated  for  convenience  as 
TZu/ctjw)  =  Zy  ( jw)  72(2  { jw) 

=  111.2(1^  +  128.5)  (jo)  +  188.5 _ 

(ju)  +  0.3397  +  jll.30)(jco  +  82  ±  jll9.8) 

(4-35) 

Using  the  straight  line  approximations  (Ref  6:255), 
the  Eq  4-35  denominator  presents  positive  slopes  of  40  dB 
while  the  numerator  gives  slopes  of  negative  80  dB  which 
results  in  a  second  order  response.  The  denominator  effects 
do  not  become  prominent  until  20  Hz.  To  explain  the 
"second  order  response"  measured  at  CGTIF  (Ref  19)  the  Z 


response  (SSP  Isolation),  Figure  4.2  is  analyzed.  A 
-30  dB  slope  approximation  is  measured  by  FJSRL  engineers 
(Ref  19)  for  the  vertical  direction.  The  two  second  order 
responses  theoretically  at  1.3  Hz  and  4.0  Hz  or  measured 
at  1.2  Hz  and  3 . 5  Hz  (Ref  19)  cause  large  magnitude  overshoots 
which  would  distort  any  straight  line  approximation  over 

such  a  narrow  bandwidth. 

Based  on  the  frequency  analysis,  the  active  LQ 
controller  synthesis  first  considers  the  second  stage 
'  dynamics.  The  reduced  order  controller  is  compared  to  a 

truth  model  representing  the  SS?  dynamics,  neglecting  pier 
dynamics.  Also  from  a  classical  controls  viewpoint,  the 
active  controller  design  is  a  difficult  problem  with  the 
small  phase  and  gain  margins  theoretically  present  for  the 
j  full  SSP  isolation  system. 
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Introduction 


V.  Active  Controller  Design 


An  active  digital  controller  is  designed  using  Linear 
Quadratic  (LQ)  synthesis  techniques  based  on  several  alternate 
descriptions  of  the  SSP  dynamics.  These  alternatives  arise 
by  considering  the  SSP  isolation  in  the  passive  response 
analysis  along  with  the  dominance  of  the  upper  level  resonances. 
The  use  of  the  upper  level  actuators  is  expanded  as  a  serious 
controls  issue  by  including  additional  active  control 
actuators  between  the  lower  level  and  the  pier. 

Each  alternative  is  developed  using  standard  state 
variable  notation  and  partitioned  submatrices  from  the 
differential  equations  written  to  describe  the  SSP  dynamics 
in  Chapter  III.  A  simulation  of  isolator  misalignments  is 
included  to  evaluate  translational  and  rotational  cross 
coupling  effects. 

A  linear  quadratic  cost  function  and  control  law  is 
defined  to  design  a  digital  controller  to  meet  the  required 
SSP  performance  specifications.  A  general  procedure  is 
outlined  for  controller  design  process  along  with  design 
policy  based  on  general  SSP  dynamics  constraints.  Each 
controller  is  then  designed  and  analyzed  with  classical 
theory  cited  to  reinforce  the  digital  performance.  Now,  the 
SSP  state  representation  is  discussed  for  the  continuous 
time  domain. 
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SSP  state  Space  Representation 


The  dynamics  of  the  SSP  describe  a  linear  time  invariant 
system  using  standard  state  variable  notation  (Ref  6:27,  12) 
as 


X(t)  =  A  X{t)  +  B  U(t) 
X(t)  =  C  X(t)  +  D  U(t) 


(5-1) 

(5-2) 


where  X ( t ) , 


state  vector  derivative  with  respect  to 
time,  column  vector 


X(t),  state  vector  column  vector 


A,  square  matrix,  order  equal  to  dynamics 
states  of  system  (n) 

B,  input  control  matrix,  n  X  r,  where  r  is 
number  of  controls,  n  is  order  of  system 

C,  output  matrix,  p  X  n,  where  p  is  number  of 
desired  outputs,  n  is  order  of  system 

D,  feed  thru  matrix,  p  X  r,  represents  controls 
appearing  in  output 


All  the  time  dependence  notation  is  later  dropped  for 
convenience  but  is  still  implied.  All  derivatives  are  with 
respect  to  time,  an  inertial  reference  frame,  with  constant 
mass,  nominal  perturbed  position  and  initial  conditions  as 
indicated.  The  SSP  state  vectors  and  the  A  system  matrix  are 
now  defined  in  partitioned  form  using  SSP  generic  state 
variables  as 
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where, 

X  ,  X  ,  and  X  are  generic  state  variables  for  translation 
^  y  ^  motion 

X  ,  X.,  and  Xg  are  generic  state  variables  for  rotation 
^  motion 


where  A,  system  matrix 

■  ,  Ay  ,  ,  represent  translation  dynamics 

A^  ,  A^  ,  Ag  ,  represent  rbtation  dynamics 

A  g  ,  A  g  ,  represent  translation-rotation 

.  ^  cross  coupling  dynamics 

Z\l)  '  Z<|) 

A,  ,  A.  ,  represent  rotation-translation 

\1)Z  <PZ  t  c  t  ■ 

cross  coupling  dynamics 

*0x  '  ^0x 

0  ,  represents  zero  matrices 
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For  each  candidate  controller  proposed,  various  A 
partitions  and  state  vectors  are  rewritten  into  either 
reduced  order  controller  models  or  partitions  representing 
a  specific  controller  loop  (really  a  single  degree  of 
freedom) .  The  state  vectors  and  A  matrix  partitions  are 
expanded  in  Appendix  J.  The  A  matrices  based  on  coeffi¬ 
cients  specified  in  Appendix  F  and  G,  are  rewritten  in  state 
form  for  Appendix  J.  To  simulate  the  cross  coupling  terms, 
isolator  misalignments  are  simulated  under  random  isolator 
placement  deviations  of  plus  or  minus  a  one-half  inch  over 
a  uniform  distribution.  Cross  coupling  coefficients  are 
also  given  in  Appendix  J  along  with  the  simulated  isolator 
placement  distances. 

As  each  controller  design  is  discussed,  the  B,  C,  and  D 
continuous  time  matrices  are  derived  for  the  LQ  synthesis. 

Root  locus  and  Nichols  plots  are  made  from  TOTAL  using 
classical  techniques  to  predict  controller  performance  and 
stability.  In  each  case  a  basic  controller  structure  is 
developed  along  with  a  truth  model  structure  reflecting 
cross  coupling  effects. 

Six  basic  controllers  are  developed,  one  for  each  of 
the  six  degrees  of  freedom  -  three  translational  (X,  Y,  Z) 
and  three  rotational  {ip ,  <p,  C).  Six  independent  control  loops 
are  derived  using  second  and  first  stage  (SSP)  isolation. 

Truth  models  are  derived  using  particular  controller  loops 
and  the  appropriate  simulated  cross  coupling  terms.  Reduced 
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order  models  are  derived  from  the  SSP  isolation  configurations 
by  recognizing  the  dominance  of  the  second  stage  isolation 
dynamics  from  the  transmissibility  and  pole-zero  analysis 
in  Chapter  IV.  Next  the  LQ  design  philosophy  is  specified 
along  with  a  description  of  the  proposed  control  law  structure. 

Linear  Quadratic  (LQ)  Synthesis 

Of  the  numerous  digital  controller  design  methodologies 
available,  linear  quadratic  synthesis  offers  an  optimum 
controller  derived  on  the  basis  of  system  state  description 
and  a  direct  interaction  with  a  performance  index.  Classical 
theory  (Ref  6:486)  describes  the  performance  index  and  cost 
weighting  matrix  applied  to  continuous  time  controllers. 

A  LQ  controller  guarantees  stability  for  any  choice  of 
weighting  matrices,  using  classical  pole  placement  techniques 
(Ref  12) .  A  digital  controller  is  proposed  for  the  SSP 
because  of  the  system  flexibility,  expandability  and  desirable 
control  characteristics  not  found  in  analog  controllers.  In 
either,  the  continuous  or  digital  controller,  the  difficult 
Riccati  (Ref  6:495)  equations  are  solved  for  the  optimum 
feedbac)c  solution.  The  digital  solution  is  usually  easier 
to  obtain  than  the  continuous  solution  (Ref  12) .  Further 
theoretical  development  on  the  LQ  synthesis  for  digital 
controllers  is  covered  in  reference  12.  The  LQ  synthesis 
used  in  this  study  is  a  proportional  integral  (PI)  controller 
which  guarantees  a  type  1  property  desired  for  disturbance 
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rejection  (Ref  12,  6:176).  Integration  of  the  controller 
with  the  SSP  physical  system  is  shown  in  Figure  5.1.  The 
LQ-PI  digital  controller  structure  is  based  on  a  control 
difference  or  "Pseudo-Rate"  (Ref  12)  whose  characteristics 
have  a  cost  weighting  matrix  assigned  to  the  system  state 
representation  which  is  augmented  to  the  control  input  rate. 
This  LQ-PI  design  assures  zero  steady  state  tracking  and 
a  non-zero  setpoint  (Ref  12) .  The  non-zero  setpoint  requires 
the  control  input  to  keep  the  SSP  in  the  nominal  position 
for  which  this  study's  dynamic  model  is  valid.  Secondly, 
the  zero  steady  state  tracking  works  to  keep  the  accelerations 
and  tilts  for  each  axis  within  the  performance  required  by 
the  specifications.  The  LQ-PI  controller  structure  with 
the  SSP  dynamics  is  shown  in  Figure  5.2A.  The  discrete 
controller  representation  is  shown  in  Figure  5.2B. 

The  LQ-PI  controller  design  is  driven  by  a  cost 
function  described  as 


J 


— T  —  — T  — 

[y  y  y  +  u  u  u 

C  C 


+  U  Uj.  U]dt 


(5-5) 


where  J , 
Y, 

^c' 

U, 


quadratic  cost  function 
system  output  desired 
output  cost  weighting  matrix 
control  input  vector 
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U^,  control  cost  weighting  matrix 

A 

U,  control  input  rate  vector 

U^,  control  rate  cost  weighting  matrix 
T 

(  denotes  transpose) 

The  quadratic  cost  function  in  Eq  5-5  provides  the 
pseudo-rate  constraint  by  the  cost  parameter.  Since 
constant  gain  matrices  are  desired  for  computation 
simplification,  the  integration  is  carried  out  ignoring 
a  final  transient  in  system  performance  and  considers  the 
time  interval  to  approach  infinity  for  the  steady  state 
system  performance. 

The  digital  Riccati  (Ref  12)  of  Eq  5-5  specifies  G* 
the  optimum  feedbaclc  gain.  This  feedback  gain  is  shown 
in  Figure  2A  separated  into  two  separate  partitions  as 

G*  =  [G*  1  G*  ]  (5-6) 

C  Cl  C  2 


The  actual  controller  gain  and  of  Figure  5. 


2B 


is  related  to  G*  by 
c 


IG*  ] 


=  [G  G  ) 

Cl  C2  ' 


(5-7a) 


‘"x  “'z  I 


>-l  B, 


(5-7b) 


[Kjj  ]  l~] 


(5-70 
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where 


'Z 

I 

B. 


'D 


D 


D 


state  feedback  matrix 
output  feedback  matrix 
state  transistion  matrix 
identity  matrix 
discrete  input  matrix 
discrete  output  matrix 
discrete  feedforward  matrix 


U,  X,  <5,  have  discrete  time  argument  t^  implied.  These 
various  forms  are  presented  to  simplify  the 

I 

discussion  later  in  the  controller  design  process. 

The  feedback  gains  in  and  would  be  actually 

implemented  in  a  real  time  SSP  controller,  while  II  in 

Eg  5-7c  represents  the  digital  approximation  of  a  continuous 

time  system  for  a  given  sampling  frequency.  The  T  symbol 

in  Figures  5.2A  and  5.2B  is  the  time  delay  between  sample 

periods  and  in  a  digital  controller  implementation  is  a 

value  held  (delayed)  in  computer  memory.  Such  diagrams  can 

be  confusing,  because,  for  instance,  a  value  for  the  difference 

of  X(t^)  and  the  output  of  a  delay  is  indicated  for  Figure 

5.2B.  This  difference  is  equal  to  the  value  of  X(tj^)  at 

the  time  t.  minus  the  value  of  X(t.  )  -  the  output  of  the 
1  1"  1 

delay  from  the  last  sample  time. 

The  digital  controller  structure  in  Figure  5.2B 
represents  the  following  control  input  expression: 
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U  (t.)  =  u  (t.  )  -  Ky[X(t.)  -  X(t.  )] 

Ultj  )  (5-8) 

1“  1 


The  closed  loop  eigenvalues  of  the  system  shown  in 

Figure  5.2B  are  ($  -  B  G*  )  which  is  shown  from  the  following 

u  o 

discrete  state  difference  equations. 


+  Bo(t.)[GJ(t.)  X(t,)l 


(5-9a) 


(5-9b) 


If  the  system  described  by  the  state  transition  matrix 
is  stable,  then  the  closed  loop  eigenvalues  of  the  LQ-PI 
structure  are  stable  (Ref  12) .  Stability  is  addressed  later 
by  controllability  and  observability  issues  in  a  general 
design  approach. 


Digital  Controller  Design  Assumptions 

Several  assumptions  in  the  controller  design  are  made 
in  the  SSP  configuration,  controller  performance  evaluation, 
and  digital  constraints.  First  the  SSP  configuration  is 
discussed . 
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A  full  representation  of  the  SSP  would  include  upper 
and  lower  level,  and  seismic  block  (pier)  structural  resonances, 
as  well  as  the  first  and  second  stage  pneumatic  isolator 
and  foundation  soil  dynamics.  By  assumption,  the  structural 
resonant  modes  of  the  SSP  are  neglected  in  this  study. 
Similarly,  the  pier  dynamics  by  actual  transmissibility 
measurements  appear  to  be  above  the  20  Hz  upper  active 
controller  bandwidth.  Theoretical  pier  investigations  in 
Chapter  IV  do  place  pier  isolation  effects  at  least  on  the 
upper  edge  of  controller  bandwidth  (Ref  1) .  Certainly 
seismic  disturbances  would  also  need  to  be  modeled  if  a 
full  truth  model  was  intended.  A  truth  model  is  a  mathematical 
representation  which  attempts  to  describe  the  known  system 
characteristics  of  a  dynamic  system.  This  study  is  concerned 
with  the  goal  of  extending  the  -40  dB/decade  (Ref  14:5) 
transmissibility  to  lower  frequency  bounds  than  presently 
possible  with  passive  isolation.  This  study  considers  the 
SSP  truth  model  as  the  dynamics  of  the  first  and  second 
stage  pneumatic  isolator  dynamics  with  cross  coupling  effects 
generated  by  possible  simulated  random  isolator  misalignments. 

Neglecting  pier  dynamics  poses  two  significant 
analytical  dilemmas  for  the  controller  design.  First,  if 
the  pier  adds  no  dynamics,  it  becomes  part  of  the  Earth 
reference  frame,  or  an  inertial  frame.  The  actuator  must 
be  connected  to  some  structure,  the  pier  in  this  case. 

Secondly,  the  actuators  being  connected  to  the  upper  level, 
allow, in  addition, a  direct  feedthrough  force  appearing 


directly  as  an  acceleration  (specific  force)  in  the 
translational  motion  and  as  a  torque  in  the  rotational 
motion  performance  measurements.  These  points  are 
brought  out  in  each  controller  design  and  assuming  no  pier 
dynamics  gives  an  optimistic  controller  environment.  If 
pier  dynamics  are  allowed,  they  would  appear  transformed 
directly  on  the  upper  level.  A  transformation  would 
occur  simply  because  the  actuator  force  is  now  with  respect 
to  the  pier  coordinate  frame,  a  non-inertial  reference. 

All  expressions  derived  in  Chapter  III,  Appendix  F  and  G 
are  specifically  derived  for  the  equations  of  motion 
describing  the  six  degrees  of  freedom  of  the  upper,  lower 
level  and  pier  centers  of  mass  with  respect  to  the  Earth 
as  an  inertial  reference  frame.  For  the  digital  controller 
design,  the  state  representation  of  the  SSP,  considers  the 
pier  as  an  inertial  reference  for  the  actuators,  and  includes 
only  the  second  and  first  stage  isolation  dynamics. 

Controller  actuators  are  proposed  between  the  upper  level 
and  pier  (Ref  14) .  Lower  level  actuators  are  added  for  this 
study  based  on  physical  intuition  that  more  actuators 
could  better  "steady"  the  lower  level  and  the  SSP. 

To  use  the  LQ  synthesis  technique,  alternate  represen¬ 
tations  of  the  SSP  dynamics  are  developed  for  candidate 
digital  controllers  by  describing  each  in  the  state 
notation  of  Eqs  5-1  and  5-2.  Cost  analysis  is  also  proposed 
for  each  in  the  forms  needed  by  Eq  5-5 .  Actual  controller 
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design  is  done  using  a  program  CGTPIF  (Ref  7)  which  has 
an  option  which  accepts  continuous  time  system  matrices 
(A,B,C,D) .cost  criteria  and  yields  a  digital  controller 
as  shown  in  Figure  5.2B.  Figure  5.2B  is  a  discrete 
representation  of  the  SSP  dynamics  and  control  law.  No 
digital  approximations  are  made  for  the  dynamics  of  the 
analog  to  digital  (A/D)  converters.  The  sampling  frequency 
is  200  Hz  which  is  selected  to  be  ten  times  20  Hz,  the 
highest  frequency  limit  on  the  controller  bandwidth. 

The  CGTPIF  digital  controller  design  software  has 
options  intended  for  time  domain  analysis  and  is  adapted  for 
LQ-PI  controller  feasibility.  No  frequency  domain  analysis 
with  a  Bode  plot  or  PSD  options  is  available.  To  evaluate 
the  LQ-PI  feasibility,  passive  isolation  time  responses 
are  made  using  TOTAL  for  a  unit  step  input  and  compared  to 
CGTPIF  time  responses  with  initial  conditions  appropriate 
to  the  TOTAL  analysis.  Early  efforts  on  the  FJSRL  Iso-Pad 
controller  design  (Ref  4)  used  unit  step  time  responses  over 
a  four  second  settling  time  criteria  to  meet  controller 
specifications.  Approximations  are  made  for  the  resulting 
exponentially  damped  sinewave  unit  step  response.  First, 
a  unit  step  can  be  scaled,  as  was  done  in  Chapter  IV 
transmissibilities  studies,  as  long  as  the  system  operation 
remains  linear.  The  unit  step  could  be  considered  as  an 
one  g  input  disturbance  for  translation  control  and  as  one 
arcsecond  tilt  disturbance  for  rotation  control  analysis. 
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These  inputs  can  be  scaled  appropriately  once  true 
disturbance  magnitudes  are  known  without  reaccomplishing 
the  controller  simulations.  Secondly,  the  damped  frequency 
is  an  approximation  of  the  acceleration  output  frequency 
content  for  the  translation  controllers.  The  tilt  specifica¬ 
tions  are  given  in  time  response  constraints  and  are 
evaluated  directly  from  the  CGTPIF  analysis  for  the  rotation 
controllers.  The  controller  design  approach  using  CGTPIF 
is  now  discussed. 


Controller  Design  Approach 

The  SSP  digital  controller  is  approached  as  six 
individual  controllers  -  one  for  each  degree  of  freedom. 
Actually,  only  four  loops  are  needed  because  the  X  and  Y, 
and  ip  and  (p  loop  pairs  are  identical. 

Within  each  loop,  two  system  A  matrices  are  possible 
by  considering  the  full  SSP  isolation  dynamics  or  a  reduced 
order  model  using  only  the  upper  level  dynamics.  In  the 
SSP  isolation  controller  structure,  three  further  permutations 
are  possible  using  combinations  of  the  upper  and  lower 
level  actuators.  These  permutations  result  in  one  multiple- 
input  multiple-output  (MIMO)  system  or  two  single-input 
single-output  (SISO)  systems  -  the  upper  or  lower  level 
actuators  for  control  inputs.  The  translation  (X,  Y,  Z) 
controller  has  a  direct  feedthrough  of  the  actuator  (specific 
force  or  acceleration)  along  with  state  combinations  to  the 


output  (upper  level)  acceleration,  while  the  rotation 
controller  (ijj,  (p,  G)  output  tilt  is  a  system  state  with 
no  feedthrough  contribution. 

In  summary,  four  configurations  are  analyzed  for  the 
two  basic  controllers  in  each  direction  designated  as 

1.  SSP  Isolation  (MIMO)  -  Upper  and  Lower  Level 

Actuators 

2.  SSP  Isolation  (SISO)  -  Upper  Level  Actuators 

3-  SSP  Isolation  (SISO)  -  Lower  Level  Actuators 

4.  Second  Stage  Isolation  (SISO)  -  Upper  Level 
Actuators 

Controller  costs  are  specified  for  Eq  5-5  in  each 
controller  design.  The  actuator  cost  weight  U  is  derived 
using  guidelines  (P.ef  11)  for  a  control  weighting  equal  to 
the  reciprocal  of  the  maximum  actuator  force  squared.  The 
maximum  SSP  force  available  is  ten  pounds  (Ref  14:24)  which 
yields  a  0,01  weighting  for  U^.  The  control  rate  U^  is 
assumed  to  be  0.1,  as  ten  times  the  actuator  cost  weight. 

The  output  cost  is  varied  from  1  to  100  with  the  results 
noted  in  each  controller  design  section.  The  cost  weighting 
matrices  have  zeros  added  for  the  MIMO  and  SISO  configurations 
to  account  for  input  control  configurations,  no  cross 
cost  considerations,  and  a  desired  system  output  penality. 

Controllability  and  observability  analysis  is  theoretically 
commented  on  in  each  controller  section  using  basic  analysis 
of  the  B  input  matrix  and  C  output  matrix  for  controllability 
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and  observability  measures  directly.  The  results  in  each 
case  is  verified  using  a  software  package  "MIMO"  (Ref  13) 
but  is  not  listed  in  the  study. 

All  configurations  suffer  from  being  weakly  controllable 
which  is  a  direct  function  of  the  coefficient  value  in  the 
B  input  matrices  being  either  scaled  by  large  mass  (trans¬ 
lation)  or  moment  of  inertia  (rotation) .  For  example, 

matrix  for  X  controller  is  scaled  by  the  upper  and  lower 
masses.  The  factor  l/m^  scales  Eq  F-5  by  0.0019  which  is 
evident  in  the  last  row  of  and  B^^  in  the  X  -  Controller 
(SSP  Isolation) .  To  avoid  the  controllability  problem, 
the  candidate  controller  design  uses  the  B  input  coefficient 
resulting  from  the  dynamics  equations  in  Appendix  F  and  G, 
but  the  CGTPIF  simulation  used  values  scaled  to  one. 

Physically,  any  input  actuator  force  is  attenuated  by  such 
a  mass  scaling  factor  and  is  describing  the  reality  that  a 
small  force  does  little  to  induce  an  acceleration  or  tilt 
to  a  large  mass. 

Observability  theoretically  is  the  most  serious  control 
problem  for  the  SSP.  Only  the  second  stage  isolation 
controller  types  have  outputs  using  upper  level  states. 

The  translation  controllers  do  not  have  their  acceleration 
outputs  as  distinct  system  state  while  the  tilt  output  of 
the  rotation  controllers  is.  The  SSP  isolation  type  controller 
has  no  outputs  for  the  control  law  taken  from  the  lower  level. 
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Unfortunately,  the  lower  level  is  physically  inaccessible 
and  instrument  monitoring  is  unlikely,  so  the  SSP  isolation  - 
MIMO  system  has  observability  problems  as  will  be  demonstrated. 

A  controller  is  now  designed  for  each  of  the  four 
basic  configurations  in  each  direction  (degree  of  freedom) . 

A  root  locus  and  Nichols  plot  is  included  to  give  insight 
into  the  controller  performance  based  on  classical  concepts. 

The  unit  step  responses  from  the  TOTAL  and  CGTPIF  simulation 
programs  are  included  for  performance  evaluation.  Digital 
controller  designs  for  successful  (with  Riccati  solutions) 
CGTPIF  simulations  are  listed  in  Appendix  K  by  direction  loop 
name  conf iguration  type . 
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X  -  Controller  (SSP  Isolation) 

The  truth  model  for  the  X  controller  is  rewritten 
from  the  partitions  of  Eqs  J-8A  and  J-8F  for  the  X  state 
variable  and  cross  coupling  terms  as 


*1 

" 

u 

x 

"0 

(5-10) 

(5-11) 


The  X  state  space  model  using  SSP  dynamics  for  X  loop, 
with  no  rotational  cross  coupling  (no  truth  model) ,  reduces 
Eq  5-10  to 


J2. 

X. 


A  X  +  B  U 
XX  XX 


(5-12) 


Using  the  state  notation  from  Eq  J-2  and  the  A  matrix 
values  from  Eq  J-10,  the  X  state  representation  is  written 
completely  as 


•  ' 

• 

o 

Xi 

0.0000 

1.0000 

0.0000 

0.0000 

o 

X2 

-2740.0000 

-6.7424 

1140.6054 

1.7230 

o 

Xa 

0.0000 

0.0000 

0.0000 

1.0000 

o 

X4 

631.6547 

1.5080 

-631.6547 

-1.5080 

L 

Xi 

X2 

Xa 

X4 


(5-13) 
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The  root  locus  of  the  X  open  loop  is  obtained  by 
considering  the  X  transmissibility  for  the  SSP  isolation 
as  a  unity  feedback  continuous  time  system.  The  root 
locus  is  shown  in  Figure  5.3.  Four  complex  poles  are 
located  near  the  imaginary  axis  with  two  poles  in  the  left 
half  plane  (LHP)  (Ref  6) .  The  lower  frequency  complex 
poles  pair  seek  out  the  two  system  zeros  while  the  higher 
frequency  poles  cross  the  imaginary  axis  into  the  right 
half  plane  (RHP)  along  a  vertical  asymptote  parallel  to  the 
imaginary  axis.  The  cross  effect  between  the  two  pair  of 
complex  poles  is  a  stability  consideration  because  of  the 
SSP  isolation  minimum  phase  and  gain  margins.  All  root 
locus  for  the  SSP  isolation  (translation  and  rotation) 
dynamics  have  the  same  general  shape  and  analysis  comments. 
The  poles  and  zeros  are  repeated  from  Eq  4-35 r 


Poles: 

Pi.  2 

=  -0.4079  ±  jl823 

Ps  ,  It 

=  -3.717  ±  j55.00 

Zeros: 

= 

-397.9 

^2  = 

-418.9 

A  Nichols  plot  for  the  open  loop  transfer  function  is 
in  Figure  5.4  which  shows  how  the  5,29  dB  (6.5  Hz)  gain  margin 
and  5.39®  (4.5  Hz)  phase  margin  affect  controller  stability. 

A  unit  step  continuous  time  response  is  shown  in  Figure  5.5 
with  a  ±0.20  damped  sine  wave  at  3.13  Hz  (on  a  bias  of  one) 
approximated at  four  seconds  for  comparison  to  the  CGTPIF 
unit  step  responses.  Other  characteristics  of  Figure  5.5 
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X  -  RESPONSE(SSP  ISOLATION) 


Fiqure  5.3  X  -  RGSoonse  (SSP  Isolation)  -  Root  Locus 
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Figure  5-4  X  -  Response  (SSP  Isolation)  -  Nichols  Plot 


RESPONSECSSP  ISOLRTION) 
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Figure  5.5  X  -  Response  (SSP  Isolation)  -  Time  Response 


are  summarized  in  Table  4.4,  under  the  TX^^pd)  response. 

Various  SISO  type  controllers  are  tried  using  the 
following  actuator  configurations  as 


Upper  Level  Actuator: 


0.0000 

0.0000 

0.0000 

0.0019 


D, 


0.0000 

0.0019 


(5-14) 


Lower  Level  Actuator: 


®X^X 


0.0000 

0.0034 

0.0000 

0.0000 


If^] 


°x 


0.0034 

0.0000 


(5-15) 


The  output  is  the  same  in  both  configurations: 


Y2 


[631.6547  1.5080  -631.6547  -1.5080] 


Xz 


X3 

X4 


(5-16) 
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F„  and  f  are  the  U  terras  in  Eq  5-12  which  represent 
XX  X 

the  X  upper  and  lower  level  actuator  resultant  forces.  The 
lower  actuators  are  indicated  by  f  in  Figure  3,4.  Resultant 
forces  from  these  actuators  use  a  subscript  for  the  summation 
orientation  on  the  SSP.  No  solutions  are  possible  from 
CGTPIF  for  these  SISO  configurations,  even  with  scaled 
to  one  to  account  for  the  weak  controllability.  These 
controllers  are  not  fully  observable  and  controllable.  More 
design  analysis  is  done  with  SISO  configurations  for  the  Z 
controller , 

A  MIMO  configuration  is  tried  using  the  upper  and  lower 
actuators  with  the  only  output  being  the  upper  level  accel¬ 
eration,  which  is  a  combination  of  system  states  and  not 
fully  observable.  The  and  C„  matrices  are  listed  below: 


0.0000  0.0000 

0.0034  0.0000 

0.0000  0.0000 

0.0000  0.0019 


(5-17) 


Yi 

0.0000 

0.0000 

0.0000 

0.0000 

631.6547 

1.5080 

-631.6547 

-1.5080 

0.0034 

+ 

0.0000 


0.0000 

0.0019 


(5-18) 
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Figure  5.6  X  -  Controller  (SSP  Isolation) 


CGTPIF  provided  a  solution  which  is  listed  as  Case  1, 

Design  A  in  Appendix  L.  An  approximate  estimate  of  the  system 
response  can  be  made  from  Figure  5.6  as  a  peak-to-peak  response  (4 
seconds)  of  ±0.02  damped  sinewave  at  3.12  Hz  on  a  bias  of 
one.  The  X  acceleration  unit  step  response  is  simulated 
by  an  initial  condition  of  0.0016  for  state  Xi .  The  Yi 
output  is  a  zero  row  because  the  lower  level  is  not  readily 
accessible.  The  input  matrix  is  scaled  to  one  to  compensate 
for  weak  controllability.  No  increase  is  possible  in  the 
output  cost  Y^  without  destabilizing  the  system  output. 


X  -  Controller  (Second  Stage  Isolation)  -  SISO 

The  upper  level  and  the  second  stage  isolation  dynamics 
partitioned  from  Eq  5-18  yield  the  reduced  order  model  as 


o 

Xi 

0.0000 

1.0000 

Xs 

0.0000 

o 

= 

+ 

X4 

-631.6547 

-1.5080 

X4 

0.0019 

(5-19) 


Yz 


[-631.6547 


-1.5080] 


[0.0019]  Fjj. 


(5-20) 


The  reduced  order  system  root  locus  is  shown  in  Figure 
5.7.  Two  complex  poles  have  their  root  locus  return  to 
LHP  an  implied  zero  in  the  far  LHP,  not  shown  due  to  plot 
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scale.  This  root  locus  plot  is  typical  for  the  second 
stage  reduced  order  models  for  all  directions.  Essentially, 
this  reduced  order  controller  does  not  go  unstable  due  to 
higher  gains,  because  the  system  eigenvalues  are  moved 
further  into  LHP.  The  second  stage  isolation  poles  and 
zeros  for  the  open  loop  transfer  function  are: 


Poles:  pi,2  =  -0.7540  ±  j25.12 

Zeros:  Zi  =  -418.9 


A  10  degree  phase  margin  can  be  approximated  from  the 
Nichols  plot  in  Figure  5.8.  Gain  margin  is  not  a  constraint 
either  in  the  Nichols  plot  or  the  root  locus  plot. 

A  unit  step  time  plot  for  the  open  loop  continuous 
time  response  is  shown  in  Figure  5.9  which  shows  +0.075 
damped  sine  wave  at  4.17  Hz  on  a  bias  of  one.  A  summary  of 
Figure  5.9  is  listed  in  Table  4.6  for  TXy^j^{t)  response. 

Case  2,  Design  A  is  the  CGTPIF  digital  controller  for 
the  X  second  stage  dynamics  reduced  order  controller.  The 
CGTPIF  unit  step  response  shows  in  Figure  5.10  (at  4  seconds) 
a  +0.08  damped  sine  wave  at  4.17  Hz  for  an  initial  condition  of 
X3  equal  to  -0.0016.  The  LQ-PI  controller  removed  the  bias 
present  before  in  Figure  5.9  produced  by  TOTAL.  For  CGTPIF, 

Bjj  is  scaled  to  one  as  weak  controllability  compensation 
and  the  output  cost  is  kept  to  one  to  prevent  instability. 
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Figure  5.8  X-  Response  (Second  Stage  Isolation)  -  Nichols  Plot 
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Figure  5.9  X  -  Response  (Second  Stage  Isolation)  -  Time  Response 


Figure  5.10  X  -  Controller  (Second  Stage  Isolation) 


Z  -  Controller  (SSP  Isolation) 


The  truth  model  for  the  Z  controller  is  rewritten  from 
the  partitions  of  Eqs  J-8c,  J-8d  and  J-8e  for  the  Z  state 
variables  and  cross  coupling  terms  as 


"  • 

•  • 

X 

BOO 

u 

z 

z 

z 

X, 

+ 

0  B,  0 

u, 

'1^ 

'p 

^  0  B , 

<t> 

(5-21) 

(5-22) 


The  Z 
reduces  Eq 


state  space  model  using  SSP  dynamics  for  Z  loop 
5-21  to 


A 


A  X 


z  z 


+ 


(5-23) 


Using  the  state  notation  from  Eq  J-4  and  A  matrix 
values  from  Eq  J-18,  Eq  5-23  is  written  as 


o 

Xg 

0.0000 

1.0000 

o 

Xio 

-589.8553 

-3.1292 

o 

Xii 

0.0000 

0.0000 

o 

Xl2 

127.9101 

0.6786 

* 

* 

0,0000  0.0000 

Xg 

230.9726  1.2254 

Xio 

0.0000  1.0000 

Xn 

127.9101  -0.6786 

Xl2 

The  Z  root  locus  plotted  in  Figure  5.11  is  analyzed 
and  commented  on  similarly  to  the  X  -  SSP  Isolation  open 
loop  root  locus. 

Poles  and  zeros  are  repeated  from  Eq  4-33: 


Poles : 

Pi  »  2 

=  -0.1883  ±  j8.422 

Pa  I  4 

=  -1.716  ±  j25.37 

Zeros: 

Z]  ,  2 

=  -188.5 

A  Nichols  plot  drawn  in  Figure  5.12  depicts  the  5.25  dB 
(3.0  Hz)  gain  margin  and  5.71°  (2.0  Hz).  A  unit  step  contin¬ 
uous  time  response  is  shown  in  Figure  5.13  with  a  ±0.60 
damped  sine  wave  at  1.32  Hz  (on  a  bias  of  one)  approximated 
at  the  four  second  time  criteria  for  comparison  to  CGTPIF. 

A  further  summary  of  Figure  5.13  is  made  in  Table  4.4  under 
the  TZy^p(t)  response. 

Several  SISO  controllers  are  proposed  using  the  following 
actuator  v-onf igurations  as. 


Upper  Level  Actuators: 


®Z^Z 


0.0000 

0.0000 

0.0000 

0.0019 


0.0000 

0.0019 


(5-25) 
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Figure  5.12  Z  -  Response  (SSP  Isolation)  -  Nichols  Plot 


RESPONSE(SSP  ISOLRTION) 


Response  (SSP  Isolation)  -  Time  Response 


Lower  Level  Actuators: 


0.0000 
0.0034 
0.0000 
0.0000 

(5-26) 

The  output  expression  is  the  same  in  both  designs: 

[Ye]  =  [127.9101  0.6786  -127.9101  -0.6786]  f Xg 

Xio 

Xii 
Xl2 

(5-27) 

Only  the  upper  actuator  configuration  (SISO)  provided 
a  digital  controller  design  from  CGTPIF.  The  output  cost 

is  varied  from  1000  to  10,000  with  the  following  performance 
results ; 

Y 

c 

1000 (Design  A) 

10 , 000 (Design  B) 

100,000 


Z  -  Acceleration 

±0.1  damped  sine  wave  (1.35  Hz,  bias  -0.5) 
±0.6  damped  sine  wave  (1.35  Hz) 
no  solution  from  CGTPIF 


0.0034 

0.0000 
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The  input  matrix  B  is  not  scaled  and  the  designs  are 

Lx 

listed  in  Case  3,  Appendix  L.  Of  all  the  SISO  controllers 
for  any  SSP  isolation,  CGTPIF  gave  a  solution  for  the  Z 
direction.  No  step  responses  are  available  but  their  form 
is  similar  to  Figure  5.14.  The  only  explanation  is  the 
double  zero  in  the  root  locus  at  -188.5  and  is  probably  a 
theoretical  dynamic  structure  not  necessarily  possible  in 
the  actual  SSP  dynamics. 

A  MIMO  design  is  built  using  the  following  B  ,  Y 

li  c 

expressions : 


0.0000 

0.0000 

0.0034 

0.0000 

0.0000 

0.0000 

0.0000 

0.0019 

Ys 

0.0000 

0.0000 

0.0000 

0.0000 

Ye 

127.9101 

0.6786 

-127.9101 

-0.6786 

Xii 

Xl2 


0.0000  0.0000 

0.0000  0.0019 

^2 

_ 

(5-29) 
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Figure  5.14  Z  -  Controller  (SSP  Isolation) 


CGTPIF  obtained  a  digital  design  listed  in  Appendix  L 


Case  3,  Design  C.  The  unit  step  acceleration  for  an  initial 
condition  on  X9  of  0.0078  is  shown  in  Figure  5.14.  The 
performance  approximated  at  4  seconds  in  Figure  5.14  is 
±0.0729  damped  sine  wave  at  1.32  Hz.  The  input  matrix  B 

Z 

is  scaled  to  one  and  could  not  be  increased  past  one. 

The  discrete  input  matrix  is  listed  in  Case  C. 


Z  -  Controller  (Second  Stage  Isolation)  -  SISO 

The  upper  level  and  the  second  stage  isolation  dynamics 
partitioned  from  Eq  5-24  yield  the  reduced  order  model  as 


0 

Xu 

0.0000 

1.0000 

Xu 

0 

X12 

-127.9101 

-0.6786 

1 

X 

[0.0019] 


Ye 


[-127.9101 


-0.6786] 


+  [0.0019] 


X12J 


(5-30) 


(5-31) 


A  root  locus  for  the  upper  level  dynamics  is  shown  in 
Figure  5.15  and  is  typical  of  the  second  stage  root  locus 
as  discussed  for  the  X  controller.  The  second  stage  isolation 
poles  and  zeros  for  the  open  loop  transfer  function  from 
Eq  4-50  are: 

Poles:  pi,2  =  -0.3393  ±  jll.30 

Zeros:  Zi  =  -188.5 
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A  10  degree  phase  margin  can  be  interpreted  from 


Figure  5.16.  Gain  margin  is  not  a  constraint  either  in 
the  Nichols  plot  or  root  locus. 

A  unit  step  time  response  for  the  open  loop  continuous 
time  system  is  shown  in  Figure  5.17  which  shows  at  four 
seconds  a  ±0.3  damped  sine  wave  at  1.79  Hz  on  a  bias  of  one. 

A  summary  of  Figure  5.17  is  listed  in  Table  4.6  for  TZy^j^(t) 
response . 

Case  4,  Design  A  is  the  CGTPIF  digital  controller  designed 
for  the  Z  second  stage  isolation.  The  CGTPIF  unit  step 
response  shows  at  four  seconds  in  Figure  5.18  a  ±0.5021 


damped  sine  wave  at  1.92  Hz  for  an  initial  condition  of  Xu 

equal  to  -0.0078.  The  designed  controller  removed  the 

bias  present  in  the  open  loop  response.  The  digital  controller 

required  B„  scaled  to  one  and  Y  equal  to  one  for  a  solution. 

^  c 


PSI  -  Controller  (SSP  Isolation) 

The  truth  model  for  the  controller  is  rewritten  from 
the  partitions  of  Eq  J-8c  and  J-8d  for  the  state  variables 
and  cross  coupling  terms  as 


(5-32) 
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Figure  5.17  2  -  Response  (Second  Stage  Isolation)  -  Time  Response 


(5-33) 


C,  X, 


D,  U, 

tjj  i|j 


The  state  space  model  using  the  SSP  dynamics 
loop  reduces  Eq  5-32  to 


X 


'll 


A,  X,  +  B,  U, 

ill  ip 


(5-34) 


Using  the  state  notation  from  Eq  J-5  and  the  A  matrix 
values  from  Eq  J-26,  Eq  5-34  is  written  as 


o 

1 

• 

Xl3 

0.0000 

1.0000 

0.0000 

0.0000 

o 

Xl4 

-1425.4116 

-7.5620 

563.1032 

2.9870 

o 

Xu 

0.0000 

0.0000 

0.0000 

1.0000 

o 

Xl6 

281.5035 

1.4934 

-284.0217 

-1.5068 

J 

* 

. 

Xl3 

Xl4 

Xu 

Xie 


(5-35) 


The  root  locus  of  the  open  loop  is  derived  from  i|s 
transmissibility  for  the  SSP  isolation  with  a  unit  feedback 
continuous  time  system.  The  root  locus  is  drawn  in  Figure 
5.19  and  is  analyzed  for  the  same  general  form  for  the  X 
SSP  isolation  controller. 

The  poles  and  zeros  are  repeated  from  Eq  4-42: 


Poles : 

Pi  »  2 

=  -0.4201  +  jl2.58 

P3,4 

=  -4.160  ±  j39.38 

Zeros: 

z  1  ,  2 

=  -188.5 
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A  Nichols  plot  is  drawn  in  Figure  5.20  which  shows 
the  6.33  dB  (4.7  Hz)  gain  margin  and  8.3°  (3.1  Hz)  phase 
margin  for  this  controller  design.  The  continuous  time 
unit  step  response  is  plotted  in  Figure  5.21  which  at  four 
seconds  is  a  ±0.3  damped  sine  wave  at  2.08  Hz  (on  a  bias  of 
one)  to  be  compared  later  to  the  CGTPIF  discrete  time 
response.  A  summary  of  Figure  5.21  is  listed  in  Table  4.4 
under  Ti(;yyp(t)  response. 

Two  SISO  controllers  are  planned  using  the  following 
actuator  configurations  as 


Upper  Level  Actuator: 


B.T,  = 
'P 


0.0000 

0.0000 

0.0000 

0.0004 


(5-36) 


Lower  Level  Actuator: 


ip  ip 


0.0000 

0.0009 

0.0000 

0.0000 


(5-37) 
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The  output  is  the  same  in  both  configurations: 


=  [0.0000  0.0000  10.0000] 


Xi3 

Xl4 

Xl5 

Xl6 


(5-38) 


Notice  the  Ye  output  is  simply  the  ip  tilt  and  also 
=  0.  The  terms  in  Eg  5-36  and  Eg  5-37  are  scaled  by 
the  X  moments  of  inertia  for  the  upper  and  lower  levels. 

Based  on  past  experience,  B^  will  be  scaled  to  one  as  was 
done  for  translational  controllers  due  to  CGTPIF  not  being 
able  to  reach  a  solution.  The  translation  analogy  to  the  large 
mass  and  small  force  or  torgue  in  this  case  parallels 
rotational  dynamics  from  the  translational  physical  realities. 
The  vector  has  upper  level  ip  torgues  (T^)  and  lower  level 
ip  torques  (t^)  .  Appendix  F,  Eg  F-14  and  F-15  relate  the 
actuators  and  torgues. 

CGTPIF  gives  no  SISO  solution  for  the  ^  controller. 

Again,  observability  and  controllability  deficiencies  plagued 
the  controller  design  process. 

A  MIMO  controller  is  proposed  using  the  following  input 
matrix  and  output  eguation: 
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B.U,  = 
'l>  'P 


0.0000 

0.0000 

0.0009 

0.0000 

0.0000 

0.0000 

0.0000 

0.0004 

(5-39) 


Yy 

0.0000 

0.0000 

0.0000 

• 

0.0000 

Yb  _ 

0.0000 

0.0000 

1.0000 

0.0000 

Xl3 

Xl4 

Xl5 

Xl6 


(5-40) 


CGTPIF  obtained  the  digital  controller  listed  in  Case 
5,  Design  A  in  Appendix  K.  The  unit  tilt  discrete  response 
is  plotted  in  Figure  5.22.  An  approximation  is  made  of 
Figure  5.22  at  4  seconds  estimating  a  +0.1760  damped  sine 
wave  at  2.0  Hz  response.  This  response  resulted  from  an 
initial  condition  for  Xi3 equal  to  one.  The  Yy  output  is  a 
zero  row  assuming  the  lower  level  is  not  accessible.  The 
matrix  is  scaled  to  one  for  weak  controllability  compen¬ 
sation  and  is  not  set  above  one,  since  the  tilt  output 
is  unstable  for  high  costs. 


PSI  -  Controller  (Second  Stage  Isolation) 

The  upper  level  and  the  second  stage  isolation  dynamics 
partitioned  from  Eq  5-40  yield  the  reduced  order  model  as 
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Figure  5.22  PSI  -  Controller  (SSP  Isolation) 


Xi5 

o 

Xl6 


0.0000 

-284.0217 


1.0000 

-1.5080 


Xl5 

Xl6 


+ 


[0.0004) 


Yb 


[1.0000  0.0000] 


Xl5 

Xl6 


(5-41) 


(5-42) 


The  PSI  second  stage  isolation  open  loop  root  locus  is 
shown  in  Figure  5.23.  This  is  a  typical  second  stage  root 
locus  as  described  for  the  X  direction.  The  open  loop 
poles  and  zeros  from  Eq  4-56  are  repeated  as: 


Poles:  pi, 2  =  -0.7499  ±  jl6.80 

Zeros:  Zi  =  -188.5 


The  typical  second  stage  10  degree  phase  margin  is  read 
from  Figure  5.24,  PSI  -  Response  (Second  Stage  Isolation)  - 
Nichols  Plot.  Neither  the  root  locus  or  Nichols  plot  indicate 
a  gain  margin  problem. 

An  open  loop  unit  step  continuous  time  response  is 
plotted  in  Figure  5.25  which  shows  at  four  seconds  a  +0.075 
damped  sine  wave  at  2,5  Hz  on  a  bias  of  one.  A  summary  of 
Figure  5.25  is  available  in  Table  4.6  for  the  Ttj^yyj^(t) 
response . 
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Figure  5.23  PSI  -  Response  (Second  Stage  Isolation)  -  Root  Locus 


Figure  5.24  PSI  -  Response  (Second  Stage  Isolation)  -  Nichols  Plot 
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Figure  5.25  PSI  -  Response  (Second  stage  Isolation)  -  Time  Response 


Case  6,  Design  A  is  the  resulting  digital  controller 
from  CGTPIF.  The  input  matrix  is  scaled  to  one  for  weak 
controllability  and  the  output  cost  is  kept  at  one  because 
increased  costs  destabilize  the  output  tilt.  The  unit  tilt 
response  is  shown  in  Figure  5.26  which  at  four  seconds  is 
a  +0.0108  damped  sine  wave  at  2.78  Hz  for  an  Xis  initial 
condition  equal  to  one. 

Theta  -  Controller  (SSP  Isolation) 

The  truth  model  for  the  3  controller  is  rev/ritten  from 
the  partitions  of  Eqs  J-8a,  J-8b  and  J-8f  for  the  0  state 
variables  and  cross  coupling  terms  as 


The  state  space  model  using  SSP  dynamics  for  the 
loop  reduces  Eq  5-43  to 


(5-45) 
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Figure  5.26  PSI  -  Controller  (Second  Stage  Isolation) 


Using  the  state  notation  from  Eq  J-7  and  the  A  matrix 


values  from  Eq  J-36,  Eq  5-45  is  written  as 


' 

• 

- 

- 

X21 

0.0000 

1.0000 

0.0000 

0.0000 

X21 

0 

X22 

-6104.7858 

-15.0218 

2560.1284 

6.1119 

X22 

0 

X23 

0.0000 

0.0000 

0.0000 

1.0000 

X23 

0 

X24 

1279.9637 

3.0557 

-1289.5446 

-3.0786 

X2^ 

.  . 

(5-46) 

The  root  locus  is  drawn  in  Figure  5.27  with  the  usual 
SSP  isolation  form  and  general  comments  made  for  the  X  SSP 
isolation  root  locus.  The  poles  and  zeros  are  repeated 
from  Eq  4-39: 


Poles:  pi»2  =  -0.8363  ±  j26.10 

P3,4  =  -8.273  +  j81.81 

Zeros:  Zi  =  -397.9 

Z2  =  -418.9 


A  Nichols  plot  is  drawn  in  Figure  5.28  showing  the 
6.22  dB  (9.8  Hz)  gain  margin  and  7.1°  (7.0  Hz)  phase 
margin.  A  continuous  time  unit  step  response  is  shown  in 
Figure  5.29  with  a  ±0.05  damped  sine  wave  at  4.17  Hz  (on  a 
bias  of  one)  approximated  at  the  4  second  comparison  time. 
More  information  on  the  time  response  is  listed  in  Table 
4.4  under  the  function  T0y^p(t). 
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Figure  5.27  Theta  -  Response  (SSP  Isolation) 


Root  Locus 


-IBO  -135  -90  -46 

PHASE  ANGLE  (DEGREES) 


Several  SISO  controllers  are  attempted  using  the  following 
actuator  configurations: 


Upper  Level  Actuator : 

0.0000 


0.0000 

0.0000 

0.0001 


[Tg] 


(5-47) 


Lower  Level  Actuator: 

0.0000 


0.0004 

0.0000 

0.0000 


[t^] 


(5-48) 


The  output  is  the  same  in  both  configurations: 


(Yi2 1  =  [0.0000  0.0000  1.0000  0.0000] 


X21 

X22 

X23 

X24 


(5-49) 


The  Y12  output  is  the  0  tilt.  No  SISO  solutions  are 
possible  from  CGTPIF.  The  input  matrix  is  scaled  to  one, 
but  observability  and  controllability  problems  prevented  an 
optimum  controller  solution. 
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A  MIMO  controller  is  specified  using  the  upper  and  lower 
actuators  with  the  output  as  the  9  tilt  for  the  following 
system  expressions; 


0.0000 

0.0000 

0.0004 

0.0000 

0.0000 

0.0000 

0.0000 

0.0001 

(5-50) 


• 

Yii 

0.0000 

0.0000 

0.0000 

0.0000 

Yi2 

0.0000 

0.0000 

1.0000 

0.0000 

X23 

X24 


(5-51) 


The  Yii  output  is  a  zero  row  because  the  lower  level  is 
physically  inaccessible. 

CGTPIF  provided  the  solution  listed  in  Case  7,  Design  A, 
Appendix  K.  An  estimate  is  made  from  the  9  unit  tilt  discrete 
time  response  in  Figure  5.30  for  the  tilt  at  4  seconds  to  be 
±0.019  damped  sine  wave  at  4.17  Hz  for  an  initial  condition 
for  state  X21  equal  to  one.  The  Bq  input  matrix  is  scaled 
to  one  for  controllability  compensation  and  the  output 
cost  remained  at  one  for  tilt  output  stability. 
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Figure  5.30  THETA  -  Controller  (SSP  Isolation) 


Theta  -  Controller  (Second  Stage  Isolation) 


The  upper  level  and  the  second  stage  isolation  dynamics 
partitioned  from  Eq  5-46 


Xaa' 

X24 

0.0000 

-1289.5446 


1.0000 

X 

_ -J 

3.0786 

X24  1 

> 

L  J 

[0.0002]  T0 


(5-52) 


Yia  = 


[1.0000 


0.0000] 


X23 


X24 


(5-53) 


The  reduced  order  system  open  loop  root  locus  is  shown 
in  Figure  5.31  and  is  typical  of  the  second  stage  isolation 
dynamics.  Poles  and  zeros  are  repeated  from  Eq  4-59  as: 


Poles:  pii2  =  -1.530  +  j35.77 

Zeros;  Zi  =  -418.9 


About  a  12°  phase  margin  is  indicated  in  the  Nichols 
plot  drawn  in  Figure  5.32.  Gain  margin  is  not  a  constraint 
in  either  the  Nichols  plot  or  root  locus. 

A  continuous  time  response  is  plotted  in  Figure  5.33  with 
a  flat  step  response  at  the  four  second  criteria  point.  A 
summary  of  Figure  5.33  is  given  in  Table  4.6  under  the  T0y^j^(t) 
response . 
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Figure  5.31  Theta  -  Response  (Second  Stage  Isolation)  -  Root  Locus 
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Figure  5.32  Theta  -  Response  (Second  Stage  Isolation)  -  Nichols  Plot 


THETfl  -  RESPONSECSECOND  STfiGE  ISOLfiTION) 
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CGTPIF  provided  a  digital  controller  design  in  Case  8, 

Appendix  K.  Design  A  is  plotted  for  comparison  in  Figure 

5.34  at  four  seconds  is  a  ±  0.0011  damped  sine  wave  operating 

at  6.25  Hz.  is  scaled  to  one  and  Y  is  limited  to  one 

0  c 

for  output  stability.  Design  B,  with  Bg  unsealed  gives  a 
±0.01  damped  sine  wave  at  6.25  Hz  on  a  -0.01  bias,  with  no 
plot  available  -  the  form  would  be  similar  to  Figure  5.34. 

Design  B  offered  no  improvement  with  increasing  . 

Summary 

Four  control  issues  result  from  the  controller  design 
effort  -  observability,  controllability  and  SISO/MIMO 
configurations.  These  areas  are  discussed  by  reviewing 
the  control  designs  for  the  translation  and  rotation  controller. 

The  translation  controller  performance  is  approximated 
by  analyzing  the  continuous  open  loop  unit  acceleration 
steps  from  TOTAL  and  the  discrete  time  step  responses  from 
CGTPIF,  a  certain  point  in  time  -  four  seconds  elapsed  response 
time.  The  unit  step  and  output  response  can  give  a  measure 
of  transmissibility ,  namely,  if  unit  input  disturbance  is 
applied  to  a  linear  system  and  the  output  is  0.01,  the  system 
offered  -40  dB  attenuation.  In  frequency  domain  analysis, 
sinusoid  inputs  would  have  to  be  assumed  for  a  rough  approx¬ 
imation.  Since  this  is  a  feasibility  study,  such  a  rough 
estimate  can  be  made.  Also  this  approximation  gives  some 
idea  of  the  frequency  content  of  the  system  output  -  but  only 
for  the  damped  sine  wave  frequency.  The  Z  -  Controller 
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Figure  5.34  THETA  -  Controller  (Second  Stage  Isolation) 


VI.  Conclusions  and  Recommendations 


Conclusions 

Analytically,  the  pneumatic  isolators  in  a  dual  isolator 
stage  configuration  have  severe  phase  and  gain  margin 
limitations  which  can  be  seen  in  any  one  of  the  six  degrees 
of  freedom  -  in  either  the  appropriate  Bode  plots,  root 
locus  plots  or  Nichols  plots.  A  single  level  of  pneumatic 
isolation  avoids  the  minimum  phase  and  gain  margin  charac¬ 
teristics  but  sacrifices  the  much  needed  passive  transmissi- 
bility  attenuation  to  acceleration  and  tilt  disturbances. 

The  lightly  damped  response  of  the  pneumatic  isolator  yields 
high  resonant  magnitude  peaks  well  within  the  SSP  operation 
of  bandwidth  and  produces  corresponding  slow  time  responses. 

The  actuator  energy  requirements  are  very  substantial 
because  of  the  large  SSP  mass  and  moment  of  inertia  scaling 
factors  present  in  all  controller  designs  evidenced  by  weak 
controllability. 

Observability  is  a  major  control  issue.  The  translation 
control  loop  output  is  an  acceleration  which  is  a  combination 
of  system  states  making  observability  difficult.  Rotational 
controllers  have  a  tilt  output  that  is  a  distinct  system  state, 
but  the  rotation  and  the  translation  controllers  both  cannot 
observe  the  lower  level  states  at  all.  No  lower  level 
outputs  are  taken  for  the  LQ-PI  control  law  because  in  the 
actual  controller  implementation  no  instrumentation  is  provided 
for  state  estimation. 
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The  rotational  performance  could  be  reviewed,  since 
the  specification  is  a  displacement  constraint  not  a  power 
spectral  density.  The  SISO  controllers  did  provide  design 
solutions  within  the  ±0.02  tilt  specification,  but  not 
realistic  in  light  of  the  truth  model.  Also  all  the  SISO 
controller  B  inputs  had  to  be  scaled  indicating  weak 
controllability.  The  Theta  -  SSP  Isolation  Controller  -  MIMO 
met  the  tilt  specification  with  the  B  matrix  unsealed, 
again  the  truth  model  simulation  deteriorated  the  performance. 
Observability  is  a  problem  for  the  lower  level  states. 

The  tilt  output  is  a  system  state  which  improved  stability 
as  opposed  to  the  translation  acceleration  output  which 
are  system  state  combinations. 

One  important  CGTPIF  application  note  did  result.  In 
cases  where  B  is  scaled  to  one,  almost  no  tolerance  is  given 
to  increasing  the  output  cost  v/hich  causes  outputs  to  go 
unstable  and  grow  with  time.  Since  this  is  a  basic  feasibility 
study,  no  further  efforts  are  made  to  delineate  the  CGTPIF 
solution  instabilities  due  to  SSP  dynamics  (observability 
and  controllability)  or  CGTPIF  numerical  difficulties,  or 
combinations  of  both.  A  quick  review  of  Appendix  K  shows 
a  wide  dynamic  numerical  range. 

As  a  precaution,  the  digital  controlled  sample  frequency 
is  increased  to  2  KHz  with  no  performance  improvement  noted 
in  the  discrete  time  step  responses  or  controller  stabilities. 
No  tabulation  or  plots  of  this  simulation  are  made. 


The  SSP  theoretical  models  should  be  taken  seriously 
since  they  accurately  predict  vertical  acceleration  resonant 
frequencies  and  account  for  the  observed  "second  order 
responses"  resulting  from  a  combination  of  high  resonant 
magnitudes  and  second  order  attenuations  discussed  in  Chapter 
IV. 


Recommendations 

1.  Further  SSP  active  controller  efforts  should  either 
verify  this  study's  dynamics  or  develop  new  dynamic  models 
that  can  be  verified  by  site  measurement.  Only  when  such  a 
dynamics  model  exists  should  a  careful  digital  controller  be 
designed,  simulated  and  verified. 

2.  Further  SSP  measurements  could  verify  or  disprove 
the  rotational  dynamics.  Cross  correlation  between  rotational 
and  translational  measurements  would  identify  the  magnitude 

of  cross  coupling  modes. 

3.  More  elaborate  active  controllers  could  be  studied 
as  discussed  in  Reference  16,  using  this  study's  dynamics 
model  as  a  starting  point. 

4.  Engineers  could  accept  the  passive  response  and 
estimate  and  predict  axis  accelerations  and  tilts,  compensating 
the  test  inertial  instrument  evaluation  profiles  based  on 
proven  parameter  estimation  schemes  (Ref  12) . 

5.  CIGTF  could  reconsider  the  passive  response  offered 
by  the  underground  test  facilities  in  abandoned  mines.  Such 
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a  natural  isolation  system  would  not  have  the  lightly 
damped,  high  resonant  peaks  that  pneumatic  isolators  do 
exhibit.  Soil  dynamics  have  higher  damping  ratios  and  could 
be  augmented  with  a  second  stage  fluid  isolation  system, 
similar  to  those  techniques  used  in  mechanical  gyroscopes. 

6.  Future  seismic  isolation  systems  could  consider 
single  level  pneumatic  and/or  fluid  isolation.  Single  level 
pneumatic  isolation  would  avoid  phase  and  gain  margin  constraints 
and  instabilities  of  dual  pneumatic  isolation.  Fluid 
isolation  would  exhibit  a  better  time  response  by  avoiding 
lightly  damped  dynamics. 
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Dimensions  as  shown. 


Figure  A. 3  Pier  and  Seismic  Block 
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Item  1 ,  CSDL  specification,  conversion  from  meru ' s 
to  arcseconds/second . 

10  ^meru  =  10  ^meru  (0.015  arcsecond/second 

meru 

=  1.5  X  10  ^  arcseconds/second 

Item  2 ,  an  error  of  one  arcsecond  on  the  Earth's 
surface  is  approximately  100  feet.  For  an  hour,  the  one 
arcsecond  error  corresponds  to  the  following  average 
velocities  and  acceleration  (expressed  in  g's). 

=  100  feet/hour  ,  average  velocity 

V 

^  V  “  av  =  100  feet  1  hour  =  3600  seconds 

^  hour  (hour)^ 

100  feet 

=  -  g,  gravity  units, 

(3600  seconds)^  32.2  feet/second^ 

_  100  feet _ 

( 3600) ^seconds^  32.2  feet 

second^ 

%  10  ^g's 

The  SSP  must  meet  10  ®g's,  to  be  ten  times  better 
than  the  measurement  environment.  TS'S  instruments  are 
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expected  to  be  ten  times  more  accurate  than  present 
instruments,  thus  10  ’g's. 

Item  3,  measurement  bandwidth  is  derived  as  follows: 

f  =  i  where  T,  period  (seconds) 

f,  frequency  (Hz) 

f  =  — -  =  10“  ^Hz 

100  sec 

for  a  one  hundred  twenty  day  period: 

_ 1 _ 

(120  days) (24  hours/day ) ( 60min/day ) ( 60sec/min) 

9.64  X  10“®Hz 
10“®H  z 
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Appendix  C.  Direction  Cosine  Matrix  Transformation 


Angular  rotation  for  the  SSP  are  taken  in  the  order 
,  (J)y  and  0y  for  the  upper  level  about  the  ei  ,62  , 
and  ea  axis  respectively.  0^  corresponds  to  an  azimuth 
angle  and  (J)y  and  are  tilt  angles  available  possibly  from 
tiltmeters . 


• 

•  • 

A* 

A. 

cos0y  sinOjj  0 

ei 

B.2 

= 

-sinOy  cos0y  0 

As. 

ea 

as 

0  0  1 

A\ 

63 

(C-1 


a  =  (0y)e 

O  O 

0)  = 


(c-2; 

(c-3: 


0  Rotation  about  63 


Ui 

1 

0 

0 

bi 

U2 

= 

0 

costp^ 

sinijjjj 

/N 

b  2 

(C 

0 

-sinij/y 

U3 

cosily 

bs 

-ub 

0) 


O 


(C-8) 


ijj  Rotation  about  bi 


Using  the  individual  direction  cosine  matrices,  a 
transformation  from  e  to  u  is  obtained  using  intintesimal 
rotations  where 


sinG  ^  0 

COS0  !\!  1 

and  any  angular  products  (for  example  )  equal  to  zero. 

Since  the  SSP  motion  is  considered  for  perturbed  small  angular 
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excursions  about  a  nominal  position  and  first  order 
approximations;  these  approximations  are  justified. 


u  =  C^®(Gy)e 


(C-9) 


u  = 


1 

0 

0 


1 

1 

-ip 


u 


0 

1 


''u 


0 

1 

0 


0 

1 


u 


u 


1 

0 


(C-10) 


u  = 


1 

-0 


-I),, 


u 


-<p 


u 


(C-11) 


Transformation  matrices  for  the  lower  level  (L)  and 
pier  (P)  are  derived  similarly,  and  are 


(C-12) 
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1 


p 


E  = 


P 

1 


-0 


e 


{C-13) 


— ub  , 

— ba 

— ae 

(0  + 

(0  + 

0) 

o  ^ 

O 

+ 

(flybi  + 

(C-14) 

,  _ n0 

Using  appropriate  direction  cosine  matrix  m  in  the 
upper  level  unit  vectors  (u)  are, 

—11©  ®  ^  ® 

10  =  +  (})y  (COSli;j^.U2  -  SintljjjUj) 

O  ^  /S  A 

+  0y(-sin4)yUi  +  cos<t>y  (sinij^yUa  +  cosiJ^yUa)) 

(C-15) 


With  small  angle  approximations  and  first  order  terms 
retained, 


— ue 
(0 


O  A  O  A  O  A 


(C-16) 


Physically  this  means  angular  accelerations  will  not  couple 
into  each  other  under  perturbed  assumptions  about  a  nominal 
position . 
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Expressions  for  the  lower  level  and  pier  angular 
accelerations  follow  the  same  derivation  and  approximations. 


_*!  ^  O  ^  O  ^  O  ^ 

0)  =  tjjj^li  +  <t)j^l2  +  0^13 


(C-17) 


_ QQ  O  A.  O  Av  O  /\ 

(iJ  =  +  'J'pP2  +  0pP3 


(C-18) 
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Appendix  D.  Perturbation  of  a  Body  Center 

To  understand  translational  and  rotational  motions  of 

the  SSP,  an  upper  level  corner  (or  suspension  point)  in  the 

ui  or  easterly  (E)  orientation  is  analyzed.  A  nominal  position 

is  defined  as  when  the  u  and  e  unit  vectors  coincide  with 

UE  and  U  the  nominal  corner  and  center  of  mass  locations 
o  o 

respectively.  A  translational  and  rotational  perturbation 
moves  the  U  frame  to  the  orientation  shown  below. 

/V 

e 


Figure  D.l,  Perturbation  of  a  Corner 


The  following  vectors  and  points  are  defined  for  the 
Figure  D.l. 


Points  - 


U  =  x^ei  +  +  z^es 


nominal  center  of 
mass 

position  for  upper 
level 

perturbed  center 
of  mass 

position  for  upper 
level 

nominal  upper  level 
corner 

position  (easterly) 

perturbed  upper 
level  corner 
position  (easterly) 


Vectors  - 


^U/Uq 


^UE/Uq 


position  vector  for  new  position  of  center 
of  mass  with  respect  to  nominal  center  of 
mass  position. 

position  vector  new  corner  displacement 
with  respect  to  nominal  center  of  mass. 


^UE/UE 

o 


position  vector  describing  new  corner 
position  with  respect  to  nominal  center 
position . 
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position  vector  locating 
corner  from  center  of  mass 
in  either  the  perturbed  or 
nominal  orientation,  respectively. 
Both  are  equal  and  constant  for 
an  assumed  solid  body 
structure . 


Ijj  -  distance  an  upper  level  (U)  ,  easterly  (E) 
^  corner  is  from  center  of  mass.  Corner 

is  assumed  in  same  vertical  plane  as 
center  of  mass. 


^UE/U  '  ^UE/l^ 


Since  the  new  corner  describes  the  needed  motion 
expressions,  two  vector  equivalences  of  are  written 

below. 

^VE/U  ~  ^UE  /V  ^VE/UE^ 

o  o  o  o 

~  ^UE®^  '*■  ^UE®^  ^UE®^  '*■  ^UE®^ 

^UE/U  ^U/U  ^UE/U  (D-3) 

o  o 

=  x^ej  +  y^e2  +  z^ej  +  ly^uj  (D-4) 


(D-1) 

(D-2) 


Using  Eq  (C-11) ,  Appendix  C, 


(C-11) 
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Eq  (d-4)  becomes 


^^'’’\jE^®^'''^%’''  ”  \je‘\i^^^  (0-6) 

Eq  (D-2)  is  the  measurements  an  inertial  instrument 
might  make  where  Eq  (D  ~6)  relates  the  translational  and 
rotational  states.  Equating  the  two  equations  by  component 
relates  measurements  and  states. 


^E  "  ^  (D-7) 

iuE  =  %  ^ 

^E  “  ^  “  ^‘*’U  (D-9) 


For  notation  reduction  UE  replaces  ^g/uo  ’ 
vector  UE  is  now  differentiated  with  respect  to  time  in 
the  Earth  frame. 


UE  -  ^,6,  +  (y„  +  +  iSjj  -  1ueJu)6,  (D-10 
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Defining  corner  perturbation  vectors  UN,  UW  and 
US  for  the  north,  west  and  south  upper  level  corners 
respectively,  position  and  velocity ^ector^are  derived 
and  differentiated  as  was  done  for  UE  and  UE. 


UN  = 

(Xu  - 

+  (Yu 

(D-11) 

ira  = 

(Xu  - 

o  ^ 

+  Yu^a 

<ZU 

^un‘^u^®^ 

(D-12) 

UW  = 

(Xu  • 

-  W«'  * 

(V  -  1 

uw®u^®^ 

+  ^2y  +  ^UW'**U^®^ 

(D-13) 

iw  = 

o  ^ 

+  (y„  -  - 

'uw®u>^^ 

+  (Zy  ■ 

(D-14) 

us  = 

(Xu 

t  lusOylSr 

+  (Yu  ■ 

+  l^u  '  ^usV®^ 

(D-15) 

~a_ 

US  “ 

(Xy 

o  ^ 

+  yyez 

+  (Zy  - 

^us'^u^^^ 

(D-16) 

These  corner  vecrnr  and  their  derivatives  are  used 
in  deriving  the  equations  of  n-otion  for  the  SSP  upper  level . 
The  lower  level  and  pier  subsysten.s  are  described  by  sr^i  at 
equations  with  the  U  subscript  replaced  by  1  and  P 
respectively.  The  coupling  of  translation  an^  rotation 
motions  is  apparent  fro.  the  UP,  ^  -d  bS  eapressrons. 
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The  second  order  approximations  developed  in  Eqs  3-7 
and  3-10  are  used  to  derive  the  spring  and  damper  coefficients 
for  the  Z  and  X  (also  Y)  directions.  References  are  cited 
indicating  information  source  in  MAC  SSP  report  (Ref  14). 


Second  Stage  Isolation 


f^  =  1.8  H  2,  vertical  natural  frequency  (Refl4;C-53) 

=  530.22  lb-sec^- Pt,  upper  mass  (Ref  14:E-33) 

k_„  =  16955.12  Ib/ft 

sz 

c„„  =  89.95  Ib-sec/ft 


=  4.0  Hz,  horizontal  natural  frequency  (Refl4:C-53) 

ksx  =  kgy  =  83728.99  Ib/ft 
^SX  ”  ^SY  ~  199.89  Ib-sec/ft 

Load  for  each  isolator  is  considered  as  one-fourth 
the  upper  level.  Damping  coefficient  was  0.03  for  vertical 
and  horizontal  directions  ( Ref 14 :C-10) . 
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First  Stage  Isolation 


=  1.8  Hz,  vertical  natural  frequency  (Refl4:C-53) 

Load  =  +  itij^  =  823.85  lb-sec^- ft  (  Ref  14:E-106) 

kp2  =  26344.68  Ib/ft 

which  exceeds  vertical  stiffness  specif ication  of  16920  Ib/ft 
(1410  Ib/in) 


Cp2  =  139.76  Ib-sec/ft 

=  3.8  Hz,  horizontal  natural  frequency 

kpx  =  kpY  =  117,412.70  Ib/ft 

which  exceeds  horizontal  specification  of  75,360  Ib/ft 
(6280  Ib/in)  (Ref  14;K-27) 


Cpjj  =  CpY  =  295.05  Ib-sec/ft 


Load  for  each  isolator  is  considered  as  one-fourth  the 
combined  mass  of  upper  and  lower  levels.  Damping  coefficient 
was  0.03  for  vertical  and  horizontal  directions  (Ref  14C-10) . 
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SSP  Moments  of  Inertia 


The  moments  of  inertia  are  listed  below,  (Refl4:D-19) 


for 

=  5095.8  ft-lb-sec^ 

~  “  2340  ft-lb/sec^ 

for  m 

Xj 

1  =  2547.7  ft-lb-sec^ 

XjZ 

^LX  ^LY  ""  1169.8  ft-lb/sec^ 

Pier  Lumped  Parameters 

All  parameters  are  simply  taken  from  MAC  study, 

( Ref  14;F-26) ,  using  damping  coefficients  resulting  from 
soil  analysis  (Reichert  Model). 

1^2;  =  1-2  X  10®  Ib/ft 

c„  „  =  9.34  X  10®  Ib-sec/ft 

G  It 

X  =  =  1-08  X  10®  Ib/ft 

c^  ^  =  c„^  =  5.41  X  10®  Ib-sec/ft 

j  A  G  X 
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These  parameters  represent  the  combined  spring, 
damper  dynamics  in  a  given  direction  for  the  total  load 
of  seismic  pier  and  the  SSP.  In  terms  of  SSP  Lumped 
Parameters  model  the  K  and  C  values  must  be  divided  by 
four  or 


GZ 

c 


GZ 


=  k2,/4  =  3.00  X  10^  Ib/ft 
=  c  /4  =  2.34  X  10^  Ib-sec/ft 


GX 


=  =  ky/4  =  2.70  X  10 7  Ib/ft 


'GX 


=  c^Y  =  ^  ^  Ib/ft 


Pier  Moments  of  Inertia 

The  pier  moments  of  inertia  are  calculated  considering 
the  seismic  mass  as  a  homogeneous  rectangular  concrete  block 
of  mass  5403.73  Ib-sec^-ft  (  Ref  14  ^’-26 )  ,  standard  moment 
of  inertia  formulas  (Ref  10:5.25).  Dimensions  of  seismic 
block  are  taken  from  Appendix  B. 


II 

H 

107,925.50 

Ib-f  t  -sec 

^PX 

234,575.42 

Ib-f t-sec 

It 

>4 

H 

329,289.80 

Ib-f t-sec 
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I  is  used  in  all  transmissibility  calculations 
since  a  higher  value  of  moment  of  inertia  represents  a 
lower  natural  frequency  and  greater  concern  to  controller 
development.  Ip^  would  concern  frequencies  outside  the 
active  controller  bandwidth,  approximately  a  horizontal 
natural  frequency  of  30  Hz. 
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Appendix  F.  SSP  Three  Degree  of  Freedom  -  Translation 

Eqs  3-21,  3-22,  and  3-23  result  from  Newton's 
Second  Law,  F  =  ma .  This  appendix  reduces  their  vector 
differential  equations  using  Eqs  D-6,  D-10,  D-11  thru 
D-16  for  the  upper  level,  lower  level  and  pier  displace¬ 
ments  and  velocities  in  the  force  equations  for  the 
e  frame - 

First,  the  vector  terms  and  common  coefficients  are 
made  in  Eqs  3-9,  3-10,  and  3-11. 
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Eq  F-1  expanded  using  transformation  from  Appendix  E. 
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^LE^L^ 

^sx 

+ 

- 

^LW®L> 

+ 

’"sz 

LS®L^ 

- 

(F-4) 


Similar  expansions  for  the  lower  level  and  pier  can  be  made 
using  transforms  from  Appendix  D .  These  expansions  are 
included  in  the  following  expressions  for  all  three  bodies 
in  the  ei  ,  62  »  and  63  degrees  of  translational  freedom. 

For  the  ei  direction  the  pure  translation  components 
and  rotation  cross  coupling  terms  have  also  been  segregated, 
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similarly  for  the  62  direction,  the  complete  translation 
equations  are  as  follows. 


00  O 


^SY^^UN 

■^’^SY^^LN  “ 

~*^SY^^UE  "  ^Uw'®U 

-  ’^SY^^UE 

^UW^®U 

■‘■^SY^^LE  "  ^LW^®L 

■*■  ^SY^^LE 
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“^HW  ~  ^HE 

(F-8) 

00  0 
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^  *^LN  - 

+4Cfy  Jp 

+4kFY 

yp 

+  Kfy 

^^PN  "  ^PS^ 

^^SY 

•^SY 

^LW^®L  ~  ^^FY 

^SY 

)  <1le  - 

1lw>®l 
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-  1UW>®L 

+ 

'‘fy 

<^LE 

-  w® 

/N 

the  63 

direction 

follow 

the 

same 

scheme 

(F-10) 


used  in  the  ei  and  ea  directions.  The  complete  translation 
equations  are  listed  below. 

Oo  O 

Vu  =  '^^sz^u  ■ 


^os^'^'u  - 

^SZ 

0 

^UE^'^U  ■■ 

^SZ^^UW 

<1ln  - 

^LS^’^L 

’^SZ'^LN  - 

- 

*^SZ<^LW  ' 

+  F  +F  +F  +F 
VE  VN  *VW  VS 


(F-11) 
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"'l^l  = 


'FZ 


^  ^sz 

4  <^FZ  ’^SZ 

o 

^p 

o 

^u 

^*^FZ  ■  ^LS^’^L  ^^FZ  ’^SZ^  *^LN  "  ^LS^'^'l 

^^FZ  ^SZ^  <^LW  "  ^LE^^L  ~  ^^FZ  ■*■  ^Sz’  ^^LW  "  ^Le'‘*'l 


'^FZ 

^^PN 

IpgJ'f’p 

+ 

^FZ 

^^PN 

^ps’ 

^FZ 

^^PW 

IpE^p 

+ 

’^FZ 

^^PW 

^PE 

^sz 

^^UN 

^US^'^U 

+ 

^SZ 

^^UN 

^US 

^SZ 

^^UW 

+ 

^SZ 

^^UW 

^UE 

( F-12) 
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NOTE:  Equations  for  63  direction  do  not  have  (lyg  - 

and  similar  terms  because  the  isolator  suspension  points, 
the  corners  of  each  level,  are  in  the  horizontal  plane 
containing  body  center  of  mass.  The  SSP  was  designed  so 
weight  could  be  added  or  subtracted  to  adjusted  center  of 
mass  for  various  gyro  test  tables  and  test  items. 
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The  following  matrices  could  be  used  to  relate  upper 
level  actuators  (or  lower)  shown  in  Figure  3.4  to  resultant 
forces  for  six  degrees  of  freedom.  An  assumption  is  made 
that  one  half  the  actuator  force  is  available  for  trans¬ 
lation  or  rotation  controller  command. 

The  vertical  actuators  are  utilized  as: 


• 

—  ■ 

- 

-  T  0 

T 

VN 

L  4 

X 

1  1 

F 

F 

VS 

L  4 

Z 

..11 

F 

VW 

0  4  L 

T 

■^Y 

F 

0  i  i 

VE 

4  L 

Similarly 


the  horizontal  actuators  are  controlled  as: 


1 

2 

1 

2 

0 

0 


(F-15) 


where,  F,  indicates  upper  actuators  (f  could  easily  be 
substituted  for  lower  level) 


L,  actuator  corner  location  from  center  of  mass. 
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Appendix  G  SSP  Three  Degrees  of  Freedom  -  Rotation 

The  mechanics  of  the  rotation  equation  derivation  are 
shown  for  the  upper  level  isolator  moments  in  vector 
notation.  Essentially  the  method  is 


I  a  = 


oo 


r  X  F. 


1 


(G-1) 


where 

^UX  '  ^UY  '  ^UZ  moments  of 
inertia  about  ei  ,62  and  §3  axes 

r^  isolator  moment  arms  (UE,  UN, 

UW,  US) 

F^  isolator  forces  for  each  corner 

sum  all  force  moment  arms  from 
UE,  UW,  UN  to  US 

The  process  proceeds  using  linear  approximations  for 
moment  arms,  namely 

^  =  luE^^ 
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UN  = 

^UN®^ 

UW  = 

US  = 

^US®^ 

No  e  components  arise  because  isolator  suspension  points 
are  assumed  to  be  in  a  horizontal  plane  containing  body 
center  of  gravity.  This  assumption  is  reasonable  for  the 
SSP  (upper  level,  lower  level)  due  to  construction  specifi¬ 
cations  and  adjustable  center  of  mass  capability  with  lead 
shot  compensation  for  gyro  table  loading.  The  pier  clearly 
violates  this  assumption,  but  will  be  dealt  with  later. 

The  process  continues  with  basically  a  cross  product 
of  all  force  terms  (right  hand  side)  of  Eq  F-1. 

oo  oo  oo 

f^ux'^'u  ^uy'*’u  = 

-dygdi)  X  (-[Cgjlfe  -  [Kg]UE) 

+  (lyjje2)  X  (-[Cg]^  -  [Kg]UN) 

+(-ly^ei)  X  (-[Cgjim  -  lKg]UW) 

+(-lygS2)  X  (-(Cglife  "  [ K g ] US ) 
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X  ([Cgllfe  +  fKglLE) 

+  (lyj^e2)  X  ([CglLX  +  (KglLN) 

+(-ly^ei)  X  ([Cg]Lw  +  (KglLW) 

+(-lyge2)  X  ([CglLS  +  [Kg]LS) 

+  T  (G-3) 

The  vector  T  is  a  torque  produced  by  the  active 
controller  actuators  and  is  specified 


T  = 


riF  +1F  +1F  +1F  )lei 

^■^ue'^ve  -^un^vn  -^uw^vw  -^us  VS'^  ^ 


+  [  (lyg)  (-Fjjg)  +  (-ly^)  (-F^)  ]e2 


UN' '  UN' 


(G-4) 


Minus  signs  are  retained  because  they  emphasize 
the  actuator  orientation.  Note  also,  the  sign  on  the 
actuators  forces  and  their  parallel  to  force  column  vectors 
in  Eq  F-4  . 

The  cross  product  operation  shown  in  (G-3)  is 
continued  for  the  lower  level . 
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+  (lj^gei)  X  (-[Cp  +  CglLE  -  [Kp  +  Kg]LE) 

+  (lj^^e2)  X  (-fCp  +  Cg]LN  -  [Kp  +  Kg]LN) 

+  (-lj^^ei)  X  (-[Cp  +  Cg]LW  -  [Kp  +  KglLW) 

+  (-lj^ge2)  X  (-[Cp  +  Cg]LS  -  [Kp  +  Kg]LS) 

+  (lpgex)  X  ([Cp]#  +  [KplPE) 

+(lp^e2)  X  ([Cp]ra  +  [KplPN) 

+  (-lpj.^ei)  X  ([Cp]#  +  [Kp]PW) 

+  (-lpge2)  X  ([CpI#  +  [Kp]PS) 


+(lypei)  X  ([Cg]UE  +  (Kg]UE) 

+  (lyj^e2)  X  ([Cg]UN  +  [Kg]UN) 

+  (-lyj^ei)  X  ([Cg]#  +  [Kg]UW) 

+(-lygez)  X  ([Cg]US  +  [Kg]US)  (G-5) 


The  moment  arms  used  for  pier  forces,  namely  l„„ei,  etc., 

Ur# 

are  certainly  bad  assumptions  since  isolator  suspension 
points  do  not  lie  in  the  center  of  mass  plane.  For  this  study 
the  approximation  is  justified  in  the  passive  response 
considerations . 
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Continuing  the  cross  product  for  the  pier,  the 
equations  are 

oo  oo  oo 

^^px'^p  ^PY'*’p  ^PZ®P^- 

+  (  Ipgei)  X  (-[C^  +  Cp]#  -  [Kg  +  Kp]FE) 

+  (  lpjje2)  X  (-[Cg  +  CplTO  -  [Kg  +  KplPH) 

+(-lp^ei)  X  (-[Cg  +  Cp]PW  -  [Kg  +  Kp]FW) 

+  (-lpga2)  X  (-[Cg  +  Cp]^  -  [Kg  +  Kp]P^) 

+(  Iggei)  X  ([Cg]OT  +  [Kg]5E) 

+(  X  ([Cg]?M  +  [KglSM) 

+  (-lg^ei)  X  ([Cg]GW  +  [Kg]^) 

+  (-lgge2)  X  ([Cg]^  +  [Kg]G^) 

+  (  ^  (tCplLls  +  [KplLE) 

+  (  X  ([Cplm  +  [KplElT) 

+  (-lj^^^e3)  X  ([Cplm  +  [Kp]T:?J) 

+  (-1^^362)  X  ([CpjLk  +  [Kp]LS) 

-  T  (G-6) 

Eq  G-4  defined  T.  The  operations  indicated  in  Eqs 
G-3,  G-5  and  G-6  are  carried  out  with  the  components  tabulated 
into  ei  ,  62  and  63  components  containing  the  pure  rotation 
and  translation  cross  coupling  terms. 
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^PS^PS^'^P  ~  ^^GZ  '‘’  ^Z^  ^^PN^PN  '*'  ^PS^PS^'^P 


^px'^p 


■''^GZ  ^^PN^GN  '*’  ^PS^GS^'^G  ■‘'  ^GZ^^PN^GN  ^PS^GS^'^G 
■‘■^FZ  ^^PN^LN  ^PS^LS^’^L  '’'  ’^FZ^^PN^LN  ■*■  ^PS^LS^'^L 
"^^GZ  ^FZ^  ^^PN  ~  ^PS^^P  ”  ^^GZ  ^  ^FZ^  ^^PN  "  ^PS^  ^P 


■*■  ^GZ^^PN 

Ips'^G 

+  '‘gz 

<lpN  - 

^PS^^G 

^FZ^^PN 

ips'  h 

+  "fz 

^^PN 

^ps^"l 

^UN^VN 

1  F 

US  VS 

(G-7) 


The  ^2  rotation  equations  are  as  follows 


^UY'^U 


"‘^SZ^^UE^UE  "''  ^UW^UW^'^U  “  ’^SZ^^UE^UE  ''’  ^UW^UW^ 


•"^SZ^^UE 

LE  ^UW 

LW^'^L 

•*■  ’"sz 

*^UE  LE  ■’■  ^UW  LW^'^L 

“‘^SZ  ^^uw 

o 

^ue'^u 

^SZ 

<iuw  - 

^UE^  ^U 

■'■^SZ  ^^uw 

-  ^ue’^l 

■*■  ^sz 

<iuw  - 

^UE^^L 

■^UE^VE 

1  F 

uw  vw 
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( 


LY^L 


■^^FZ^^LE^PE  '*'  ^LW^PW^'^P  ■*■  ’^FZ^^LE^PE  '*'  ^LW^PW^'^P 
■^^SZ^^LE^UE  ^LW^PW^'^U  ’^SZ^^LE^UE  ^LW^PW^'^'U 
“^''fZ  ^SZ^  “  ^LE^^L  "  ^^FZ  ^SZ^^^LW  ~  ^LE^  ^L 


+  :  c 


FZ 


+  c 


SZ 


^^LW  - 

^LE^  ^P 

*  "rz 

^^LW 

^le 

<1lw  - 

^le’^u 

*  “sz 

^^LW 

^LE 

)  z 


)  Z 


p 

U 


oo  o 

IpY<t'p  =  -(^GZ  ^FZ^  ^^PE^E  ■*■  ^pur’rw^'^’p  "  ^\z  ’Vz^  ^^E^PE  ■*■  ^PW^PW^'^P 

■^‘^GZ^^PE^GE  ^PW^GW^'^G  ■*■  ^GZ  ^^PE^GE  ^PW^GW^  "^G 

■^°FZ^^PE^LE  ^PW^LW^^L  ^FZ  ^^PE^LE  ^PW^LW^ ‘*’l 

"^"^GZ  '^FZ^  *^PW  ~  ^PE^^P  "  ^^GZ  ’^FZ^  ^^PW  “  ^PE^  ^P 


i^ipw 

^PE^  ^G 

'’■  ’^GZ^^PW 

^PE 

:  ^^pw 

-  IpE^^L 

■'■  ^FZ^^PW 

^PE 

F„„  - 

1  F 

:  VE 

UW  VW 
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Similarly,  the  63  rotation  equations  are  as  follows 


00 

^UZ®U 


. 

- 

^SX^^UN^UN  ^US^US^ 

^SX^^UN^UN  '''  ^US^US^ 

0  - 

+ 

U 

+ 

^SY^^UE^UE  ■''  ^UW^UW^_ 

^SY^^UE^UE  '’’  ^UW^Uw\ 

‘^SX^^UN^LN  ■‘'  ^US^Ls' 

^SX^^UN  LN  ■*■  ^US  LS^ 

+ 

+ 

0 

®L  + 

+ 

^SY^^UE^LE  ‘‘'  ^UW^LW^ 

^SY^^UE^LE  ■*'  ^UW^LW^ 

^sx^^us  ^sx^^us 


+ 

‘^SX^^US 

0 

1  )  X 

UN'  L 

■*■  ^SX^^US 

- 

°SY^^UE 

^SY^^UE 

+ 

^SY^^UE 

^SY<^UE  ■ 

^UE^HE  ^UN^HN  ^UW^HW  "  ^US^HS 
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oo 

LZ®L 

<^FX  ^SX'  ^  WlN  ■*■ 

O 

Q-r  - 

^’Vx  ^SX^  ^hl^hM  ^LsW 

+ 

+ 

<^FY  ■*■  *^SY^  ‘^LE^  ^liW^w' 

L 

(kpy  ^SY^  ■*■ 

. 

'^FX^^LN^PN  ^LS^PS^ 


^FY^^LE^PE  ^LW^Pw' 

•^SX^^LN^UN  ^LS’-US' 

+ 

^SY^^LE^UE  ^LW^UW^ 


Op  + 


’^FX^^LN^PN  ^LS^Ps’ 


’'fY^^LE^PE  ^LW^PW^ 


^SX^^LN^UN  ’*’  ^LS^US^ 


'^SY^’-LE^UE  ^LW^Uw'j 


G 


U 


^^FX 

■*■  ^sx' 

^LN^'^L  '  ‘’"fX  ^SX^  ‘^LS 

°FX 

^^LS 

^LN^^P 

^FX^^LS 

“  ^LN^^P 

°SX 

^^LS 

^LN^^U 

’^SX^^LS 

"  ^LN^^U 

(CpY 

■*■  ^SY 

^  ^^LE 

lLW>yL  - 

t^FY  '^Sy'  ^^LE 

^FY 

^^LE 

^LW^^P 

■*■  ^FY^^LE 

^LW^^P 

+  ^SY<^LE  ■  ^LW’^U  •"  ’"SY^^LE  ' 
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oo 

^PZ®P 

<^GZ  ‘^FX’  ^  WpN  ^PS^PS) 

^’^GX  ^  *  WpN 

- 

+ 

(Cqy  "'’  ^fy'  ^^pe^pe  ^  ^pi’f'  ro^ 

o 

®p- 

+ 

^’^GY  ’Vy^  ^^PE^PE  ■*■  Wpw’ 

p 

^GX^^PN^GN  ■*' 

^PS^GS^ 

• 

^Gx'^PN^GN  ■*■  ^PS^GS^ 

+ 

o 

®G 

+ 

‘^gy^^pe^ge  '*' 

^PW^GW^ 

^GY^^PE^GE  '*■  ^PW^GW^ 

-* 

^FX^^PN^LN  ■*■  ^PS^LS^ 

+ 

*^FX^^PN^LN  ^PS^LS^ 

+ 

o 

®L  ^ 

^FY^^PE^LE  ■*■  ^PW^LW^ 

^FY^^PE^LE  *  ^PW^LW^_ 

-  <^GX  '  ^Pn’^P  ■  ^’"gX  ’^FX’  ‘^PS  "  ^PN^^P 

+  '"gX^^PS  ~  ^Pn’^G  ■*■  ^GX^^PS  ■  ^PN^^G 

+  Cpjjdpg  -  IpN^^L  ’*’  ^FX^^PS  ~  ^PN^^L 

-  <^GY  "  ^FY^  <^PE  “  ^PW^^P  ■  ‘’"gY  ’"fY^  ^^PE  '  ^PW^^P 

+  ^GY^^PE  ”  ^PW^^G  ■*■  ^GY^^PE  “  ^PW^^G 

+  Cpydpg  -  l-pw^^L  ’^FY^^PE  "  ^PW^^L 

•i  TT  4.  1  F  (G-9) 

■*■  ^UE^HE  ”  ^UN^HN  ”  ^UW  HW  US  HS 
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Appendix  H .  Simple  Three  Mass  System 
Transfer  Function  Polynomial  Coefficients 


A 

(s) 

As^  +  Bs^ 

+  Cs  +  D 

"g 

(  s) 

Es^  +  Fs® 

+  Gs"  +  Hs^  +  Is^  +  Js  +  K 

where  numerator  terms  are 

A 

= 

^S^F^G 

B 

= 

C 

= 

F^S  '^S^F^; 

D 

= 

and  denominator  terms  are 

E 

= 

ii^mj^mp 

F 

= 

"bmptcp  +  Cg) 

+  it^m^(Cp  +  C^)  +  mj^mpCg 

G 

= 

mpn^(ks  +  kp) 

+  mpCgCp  +  mpmJCg 
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H  =  mp(  Cgkp  +  Cpkg)  +  n^Cg(kp  +  kg) 

+n}j(Cpkg  +  Cgkp)  +  CgCgCp  +  Itl  ^^k  g  (  +  C^) 

+n^jkg(CF  +  Cg)  +  m^Cg()<^  +  k^) 

I  =  m^kg(kF  +  K)  +  mpkgkp  +  n^kpd^  +  kg) 

■"  Vg^  + 

J  =  C 
K  =  D 

All  terms  in  the  coefficients  A  thru  J  are  dimensionally 
consistant  using  basic  units  of  L,  T,  M  or  length,  time 
and  mass,  respectively. 
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Appendix  I .  Transmissibility  Polynomial  Coefficients 


The  following  coefficients  were  calculated  from 
lumped  parameter  spring  and  damper  values  given  in 
Appendix  E  using  algorithums  developed  for  the  full,  dual 
and  single  isolation  transmissibilities . 

Full  Isolation 

The  full  isolation  transfer  function  numerator  and 
denominator  polynomial  are  expressed  in  the  notation  of 
Appendix  H 


TZy^g(ja))  polynomial  coefficients  are 

A  =  223.3061 

B  =  112,875.4751 

C  =  18,750,416.3104 

D  =  1,019,400,524.282 

E  =  1.0000 

F  =  176.7548 

G  =  23,603.6706 

H  =  209,405.7952 

I  =  16,105,109.6289 

J  =  C 
K  =  D 
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TXy^gCjw)  or  TY 


(ju))  polynomial  coefficients  are 


A  = 
B  = 
C  = 
D  = 


E  = 
F  = 
G  = 
H  = 
I  = 
J  = 
K  = 


T0y/G(3u.) 


A  = 
B  = 
C  = 
D  = 


E  = 
F  = 
G  = 
H  = 


U/G 

606.8114 

616,791.8850 

200,095,712.8612 

20,192,239,747.25 

1.0000 

108.5847 

24,277.2135 

508,769.0142 

69,169,079.5072 

C 

D 


polynomial  coefficients  are 


886.4462 
901.025.9320  ' 
292,305,120.3192 
29,497,358,961.09 

1.0000 

50.3383 

14,480.0687 

379,095.3026 
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I  =  53,235,281.9575 

J  =  C 
K  =  D 

polynomial  coefficients  are 

A  =  590.0608 

B  =  298,260.4301 

C  =  49,545,813.4351 

D  =  2,693,648,362.547 

E  =  1.0000 

F  =  93.6203 

G  =  13,351.4852 

H  =  248,416.7087 

I  =  19,281,576.4410 

J  =  C 
K  =  D 

Dual  Isolation 

The  dual  isolation  transmissibilities  polynomial 
coefficients  are  listed  using  the  notation  of  Eg  4-29. 
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polynomial  coefficients  are 


TZu/p(j<^) 

a  =  1.2919 

b  =  487.0599 

C  =  45,904.7099 

d  =  1.0000 

e  =  3.8078 

f  =  719.0572 

g  =  b 
h  =  c 


polynomial  coefficients  are 


a  =  117.2864 

b  =  37,177.3628 

c  =  2,840,484.2000 


d  =  1.0000 

e  =  173.6255 

f  =  22,471.5806 

g  =  b 
h  =  c 
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TZ  /  (jw)  polynomial  coefficients  are 
U/  G 

a  =  111.2444 

b  =  35,261.6579 

C  =  2,694,091.4752 


d  =  1.0000 

e  =  164.6773 

£  =  21,313.4451 

g  =  b 

h  =  c 


TXy/p(jw) 


polynomial  coefficients  are 


a  =  6.0610 

b  =  4950.7948 

c  =  1,010,309.3674 


d  =  1.0000 

e  =  8.2503 

f  =  3377.7858 

g  =  d 
h  =  c 
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polynomial  coefficients  are 


TXj^/p(ja)) 

a  =  143.4204 

b  =  85,704.0493 

c  =  11,393,496.8012 

d.  =  1.0000 

e  =  101.7669 

f  =  20,786.5958 

g  =  b 
h  =  c 


polynomial  coefficients  are 


a  =  .143.1919 

b  =  88,565.0200 

c  =  11,973,740.2255 


d  =  1.0000 

e  =  96.6045 

f  =  19,789.7791 

^  =  b 

h  =  fc 
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polynomial  coefficients  are 


TOy/p(ja)) 

a  =  27.6590 

b  =  22,592.4016 

c  =  4,610,434.4703 

d  =  1.0000 

e  =  18.2194 

f  =  7470.1610 

g  =  b 
h  =  c 


TOL^p(ju)  polynomial  coefficients  are 


a  =  96.5325 

b  =  57,685.1711 

C  =  7,668,666.9828 


d  =  1.0000 

e  =  35.1309 

f  =  7720.9232 

g  =  b 
h  =  c 
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™u/g<3“> 


polynomial  coefficients  are 


a  =  97.3424 

b  =  60,206.8298 

c  =  8,139,793.1175 


Cl  =  1.0000 

e  =  34.9107 

f  =  7,747.9586 

g  =  b 
h  =  c 


Ti|;yyp(j(o)  polynomial  coefficients  are 

a  =  6.9904 

b  =  2635.3787 

c  =  248,380.7404 

d  =  1.0000 

e  =  9.1608 

f  =  1,733.7900 

g  =  b 

h  =  c 
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polynomial  coefficients  are 


Ti|JL/p(  jw) 

a  =  131.1338 

b  =  41,566.6994 

C  =  3,175,845.3040 

d  =  1.0000 

e  =  86.0130 

f  =  11,278.3372 

g  =  b 

K  =  c 


polynomial  coefficients  are 


a  =  125.2327 

b  =  39,695.5854 

c  =  3,032,856.2097 


6.  =  1.0000 

e  =  85.0357 

f  =  11,142.5359 

g  =  b 
h  =  c 


/» 
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single  Isolation 

The  single  isolation  transmissibilities  are  listed  in 
standard  polynomial  transfer  function  form. 


TZy/L(ja)) 


=  0.6785ja)  +  127.9100 


+  0.6785(ja))  +  127-9100 


TZj^^p(ja,) 


=  0. 678510]  +  127.900 


(jo))^  +  0.6785(ja))  +  127.9100 


TZp/G(jaj)  =  149.9780ico  +  19,269.1221 


(jo))^  +  149.9780  (jo))  +  19,269.1221 


=  1.5079jo)  +  631.6547 


(jo))^  +  1.5079  (jo>)  +  631.6547 


TXL/p(joj)  =  1.4325ju)  +  570.0683 


(joj)^  +  1.4325(ja3)  +  570.0683 


TXp/G^^oj)  =  86.8716io)  +  17,342.2099 


(jo))^  +  86.8716  (ja))  +  17,342.2099 


TGu/L^jw)  =  3.0607jo)  +  1,282.0922 


(jo))^  +  3.0607(ja))  +  1,282.0922 


T0j^^p(jo))  =  3.0120jo)  +  1,198.6119 


(joj)^  +  3.0120ja)  +  1,198.6119 


TGp/cCju))  =  31.3219io)  +  6,252.8156 


(jo))^  +  31.3219{jo))  +  6,252.8156 
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TiI^u/l  =  1.4997ju)  +  282.6902 _ 

(jw)^  +  1.4997  (ju))  +  282.6902 

Ti^L/p(jw)  =  1.5535ju)  +  292.8440 _ 

(jw)2  +  1.5535(ju))  +  292.8440 

T4ip/G(jw)  =  81.7504j()j  +  10,503.2636 _ 

(joj)^  +  81.7504  (jw)  +  10,503.2636 
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notation  and  then  in  terms  of  SSP  isolator  placement  lengths 


from  each  body  center  of  mass,  spring  and  damper  variable. 


Zero  matrixes  are  omitted,  but  understood  from  Eg 


(J-8a) 

(J-8b) 

(J-8c) 

(J-8d) 

(J-8e) 

(J-8f) 


5-4.  Each  partition  is  represented  from  Appendix  P  and  G 
differential  equations  as  shown  on  the  following  pages. 
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*D-A11!  559  AIR  FORCE  INST  OF  TECH  NRItHT-RATTERSON  AFB  ON  SCHOO— ETC  F/a  it/t 

unclassified  AFtT/OE/CE/aiO-22 


272 


1245  0.0007  -0.1245  -0.0007 


1 


/ 


275 


276 


277 


I 


278 


1279.9637  3.0557  -1289.5446  -3.0786 


Random  Isolator  Placement 

The  preceding  A  partitions  were  calculated  using  random 
isolator  placement.  The  isolator  misalignments  were  assumed 
to  be  a  uniform  distribution  having  a  mean  of  53  inches  and 
over  a  uniform  distribution  of  plus  or  minus  one-half 
inch.  Granted,  these  misalignments  do  represent  a  "worse 
case"  from  MAG's  specification  (Ref  10:15)  of  plus  or  minus 
a  one-quarter  inch  isolator  placement. 

The  placements  were  generated  by  a  GDC  6600  FORTRAN 
function  RANF  which  returns  random  numbers  between  0  and  1 . 
The  isolator  locations  were  calculated  according  to  the 
following  equations. 

1  =  RANF  +52.5 

Placement  values  produced  for  the  SSP  isolators  are 
given  as 


UE 

= 

53.0796 

^UN 

= 

53.4504 

UW 

= 

53.2860 

^US 

= 

52.7976 

LE 

= 

52.9536 

^LN 

= 

52.5060 

LW 

= 

52.7748 

^LS 

= 

52.8048 

PE 

:= 

53.1888 

^PN 

= 

52.8816 

PW 

= 

52.6320 

^PS 

= 

53.3316 
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.X  K  -  CGTPIF  Digital  Controller  Designs 


Appendix  K  contains  the  results  of  designing  possible 
digital  LQ-PI  controllers  with  CGTPIF.  Only  controllers 
proposed  in  Chapter  V  with  solutions  are  listed  on  the 
following  pages.  Controller  variables  follow  Eqs  5-6  and 
5-7  for  Figures  5.2A  and  5.2B  in  Chapter  V.  Each  controller 
type  design  is  referred  to  by  a  case  number  from  discussions 
in  Chapter  V  for  different  actuator  configurations  and  cost 
criteria  designs  indicated  by  a  letter. 


Case  1  X  -  Controller  (SSP  Isolation) 
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Case  2  X  -  Controller  (Second  Stage  Isolation) 
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Case  3,  Design  A _ Z  -  Controller  (SSP  Isolation)  -  SISO 
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Controller  (SSP  Isolation)  -  SISO 
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magn  ittide  damped 


ron^-roller  {Second  Stage  Isolation) 
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Case  5 _ PSI  -  Controller  (SSP  Isolation) 
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